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TOXICITY  ASSESSMENTS  OF  ANTIMONY,  BARIUM, 

BERYLLIUM,  AND  MANGANESE  FOR  DEVELOPMENT 
OF.ECOLOGICAL  SOIL  SCREENING  LEVELS  (ECO-SSL) 

USING  EARTHWORM  {EISENIA  FEITDA)  BENCHMARK  VALUES 

1.  INTRODUCTION 

The  U.S.  Environmenteil  Protection  Agency  (USEPA)  is  developing  Ecological 
Soil  Screening  Levels  (Eco-SSLs)  for  ecological  risk  assessment  of  contaminants  at  Superfund 
sites.  Eco-SSLs  are  soil  concentrations  of  chemicals  which,  when  not  exceeded,  will 
tiieoretically  protect  terrestrial  ecosystems  from  imacceptable  harmful  effects.  They  are  derived 
using  data  generated  from  laboratory  toxicity  tests  with  differrait  test  organisms,  wWch  represent 
the  vast  array  of  ecological  recq)tors.  Whenever  sufficient  quantity  and  quality  of  information 
existed,  Eco-SSLs  for  soil  invertebrates  were  developed  from  studies  reported  in  literature. 
However,  insufficient  information  to  generate  Eco-SSLs  for  barium  (Ba),  beryllium  (Be), 
Manganese,  (Mn),  and  antimony  (Sb)  necessitated  standardized  toxicity  testing  to  fill  the  data 
gaps. 


This  study  was  designed  to  produce  benchmark  data  for  the  development  of  Eco- 
SSLs  for  Ba,  Be,  Mn  and  Sb  for  soil  invertebrates,  and  meet  specific  criteria  (USEPA,  2000), 
including:  (1)  tests  were  conducted  in  soil  having  physico-chemical  characteristics  that  support 
relatively  high  bioavailability  of  metals;  (2)  experimental  designs  for  laboratory  studies  were 
documented  and  appropriate;  (3)  both  nominal  and  analytically  determined  concentrations  of 
chemicals  of  interest  were  reported;  (4)  tests  included  both  negative  and  positive  controls; 

(5)  chronic  or  life  cycle  tests  were  used;  (6)  appropriate  chemical  dosing  procedures  were 
reported;  (7)  concoitration-response  relationships  were  reported;  (8)  statistical  tests  used  to 
calculate  the  benchmark  and  level  of  significance  were  described;  and  (9)  the  origin  of  test 
species  were  specified  and  appropriate. 

Several  soil  invertebrate  toxicity  tests,  for  which  standardized  protocols  have  been 
developed,  can  effectively  be  used  to  assess  the  toxicity  and  to  derive  protective  benchmark 
values  for  metals  (Stephenson  et  al,  2001;  Lokke  and  Van  Gestel,  1998).  We  used  the 
Earthworm  Survival  Test  (Greene,  et  al,  1989)  and  the  Earthworm  Reproduction  Test  (ISO 
1998;  Van  Gestel  et  al,  1989)  in  ^s  study.  These  tests  were  selected  on  the  basis  of  their  ability 
to  measure  chemical  toxicity  to  ecologically  relevant  species.  The  rq)roduction  test  was  used  to 
meet  the  Eco-SSL  requirement  of  inclusion  of  at  least  one  reproductive  component  among  the 
measurement  endpoints.  The  survival  test  was  included  to  compare  lethal  and  non-lethal 
toxicity. 

Special  consideration  in  assessing  chemical  toxicity  for  Eco-SSL  development 
was  given  to  the  effects  of  aging/weathering  of  soil  contaminants  on  the  exposure  of  relevant 
ecological  receptors,  as  commonly  occurs  at  Superfund  sites.  During  chemical  aging/weathering 
in  soil,  reduction  in  tiie  exposure  to  the  chemical  may  occur  due  to  volatilization,  microbial 
degra^tion  and  immobilization,  or  other  fate  processes  (e.g.,  photodecomposition,  hydrolysis, 
and  hysteresis,  etc.).  This  can  result  in  a  dramatic  reduction  in  the  amount  of  chemical  tiiat  is 
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bioavailable,  compared  to  tests  conducted  with  freshly-amended  chemicals  or  those  tested 
following  a  short  ^uilibration  period  (e.g.,  24  h).  Standardized  methods  for  aging/weathering 
of  chemicals  in  soil  are  not  available.  We  used  the  approach  developed  to  simulate  at  least 
partially,  the  aging  and  weathering  process  that  included  exposing  soils  amended  with  chemicals 
to  periodic  alternating  wetting  and  air-drying  cycles  for  3  weeks,  in  a  green  house. 


2.  MATERIALS  AND  METHODS 

2.1  Soil  Collection  and  Characterization. 

The  soil  used  in  diese  studies  was  Sassafras  sandy  loam  [Fine-loamy,  siliceous, 
mesic  Typic  Hapludult]  (SSL)  collected  from  a  grassy  field  (M-Field)  at  Aberdeen  Proving 
Ground,  MD.  Vegetation  and  file  organic  matter  horizon  were  removed  and  the  top  six  inches 
of  the  A  horizon  were  then  collected.  The  soil  was  sieved  through  a  5mm^  mesh  screen,  air- 
dried  for  at  least  72  h  and  mixed  periodically  to  ensure  uniform  drying,  passed  through  a  2-mm 
sieve,  thra  stored  at  room  temperature  before  use  in  testing.  Soil  was  dien  analyzed  for  physical 
and  chemical  characteristics  by  the  Cooperative  Extension  Service,  University  of  Maryland  Soil 
Testing  Laboratory,  College  Park,  MD.  Results  of  these  analyses  are  presented  in  Table  1. 


Table  1.  Physical  and  chemical  characteristics  of  Sassafras  sandy  loam  soil. 


Soil  Parameter 

Sassafras  Sandy  Loam 

Sand% 

71 

Silt  % 

18 

Clay% 

11 

Texture  % 

Sandy  loam 

K.9I  SLWTTiTtl  1 

4.27 

Organic  matter  % 

1.2 

_pH 

5.0 

2.2  Earthworm  Culture. 

Earthworms  (E.fetida)  were  bred  in  plastic  containers  filled  with  approximately 
14  kg  of  a  1 : 1  mixture  of  sphagnum  PRO-GRO  peat  moss  (Gulf  Island  Peat  Moss  Co.,  PEI, 
Canada)  and  BACCTO®  potting  soil  (Michigan  Peat  Co.,  Houston,  TX,  USA).  The  pH  was 
adjusted  to  6.26  ±  0.07  by  adding  calcium  carbonate  (pulverized  lime).  The  culture  was  kept 
moist  at  21±2°C  with  continuous  light.  Earfiiworm  colonies  were  fed  biweekly  with  dehydrated 
alfelfa  pellets  (27%  fiber,  17%  protein,  1.5%  fat;  OB  of  PA,  York,  PA)  that  were  fermented, 
dried,  and  ground  to  a  course  powder.  Cultures  were  synchronized  so  that  all  worms  used  in 
each  test  were  approximately  the  same  age.  Adult  worms,  0.3g  to  0.6g,  with  fully  developed 
clitella  were  used  for  testing. 
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2.3 


Chemicals  and  Reagents. 


Assessments  were  performed  using  sulfate  salts,  including  BaS04  (CAS  #7727- 
43-7, 97%;  stock  #13989;  lot  #110120,  Alfa  Aesar),  BeS04*4H20  (CAS  #7787-56-6, 99.99%; 
stock  #16104;  lot  #H09J07,  Alfa  Aesar),  MnS04*H20  (CAS  #10034-96-5,  ACS,  98.0-101.0%, 
stock  #33341;  lot  #118129,  Alfa  Aesar),  and  Sb2(S04)3  (CAS  #7446-32-4, 97%,  stock  #B492;  lot 
#L21I28,  Alfa  Aesar).  Additional  tests  were  done  for  Ba  and  Sb  to  determine  how  carrier  salts 
and  their  solubilities  affect  the  toxicity  to  E.fetida.  For  Ba,  these  compounds  included  BaO 
(CAS  #1304-28-5, 97%,  lot  #12101BI,  Aldrich  Chemical  Company),  Ba(N03)2  (CAS  #10022- 
31-8,  ACS,  lot  #000420,  Fisher  Scientific  Co.),  and  Ba(C2H302)2  (CAS  #543-80-6,  ACS,  lot 
#995963,  Fisher  Scientific  Co.).  For  Sb,  we  used  Sb  D-tartrate  Sb2(C4H406)3*6H20 
(CAS  #  126506-93-2,  lot  #1 1 1004-2,  Pfaltz  &  Bauer).  The  positive  control  used  in  this  study 
was  4-Nitiophenol  (CAS  #100-02-7, 98%,  lot  #6623HE,  Aldrich).  The  main  carrier  salt  control 
was  sulfate  as  CaS04*2H20  (CAS  #10101-41-4,  ACS,  lot  #C07704,  J.T.  Baker).  All  reagents 
used  in  extraction  of  chemicals  fi»m  soils  were  either  reagent  or  trace  metal  grade  and  purchased 
from  commercial  suppliers.  Purified  water  (ASTM  type  I;  American  Society  of  Tes^g  and 
Materials,  http://www.astm.oig)  obtained  using  Milli-RO®  10  Plus  followed  by  Milli-Q  PF  Plus 
systems  (NCllipore®,  Bedford,  MA)  was  used  throughout  the  studies.  Glassware  was  wash^ 
with  phosphate-free  detergent  followed  by  rinses  with  tap  water,  Milli-RO®  water,  nitric  acid  1% 
(v/v)  and  finally  with  Milli-Q®  water. 

2.4  Soil  Amendment  Procedure. 

Treatment  concentrations  for  toxicity  tests  with  all  sulfate  salts  and  Ba  oxide  were 
prepared  by  adding  test  chemicals  directly  to  SSL  soil  in  appropriate  proportions  to  achieve 
nominal  target  concentrations.  Soil  was  mixed  for  3  h  on  a  three  dimensional  rotary  mixa*. 

After  mixing,  soil  was  hydrated  with  purified  water  to  1 00%  of  the  soil  water  holding  capacity 
(WHC;  1 8%  water,  on  the  basis  of  the  dry  soil  mass)  for  toxicity  testing,  or  60%  of  the  WHC  for 
the  aging/weatheTing  procedure  (see  Paragraph  2.5).  Soil  prepared  for  toxicity  tests  was  allowed 
to  equilibrate  for  24  h  before  exposing  earthworms.  The  exception  was  soil  amended  wilh  Ba 
acetate,  which  was  incubated  for  5  days  to  allow  for  acetate  degradation  by  soil  microbes. 
Treatment  concentrations  of  Ba(C2H302)2»  Bap403)2  and  Sb2(C4H406)3  were  prepar^  by 
dissolving  appropriate  amounts  of  chemical  in  purified  water,  then  hydrating  pre-weighted 
amounts  of  SSL  soil  to  achieve  target  treatment  concentrations  in  soil  for  each  chemical, 
respectively,  at  the  required  moisture  level. 

2.5  Aging/Weatfaering  of  Amended  Soil. 

All  soil  treatment  concentrations  and  negative  controls  were  subjected  to 
simulated  aging/weathering  procedure,  which  included  alternating  wetting/air-drying  cycles  for 
3  weeks  prior  to  commencement  of  definitive  tests.  Aging/weathering  of  test  soils  was 
conducted  in  open  plastic  bags  in  the  green  house.  Soil  treatments  were  initially  hydrated  to 
60%  of  water  holding  capacity  (WHC),  and  then  allowed  to  begin  drying.  All  soil  treatments 
were  weighed  and  adjusted  to  60%  of  WHC  twice  each  week,  and  afterward  brought  to  100%  of 
WHC  (18%  water,  on  the  basis  of  the  dry  soil  mass)  for  initiation  of  bioassays.  A  separate  study 
was  conducted  using  Mn  as  a  model  chemical  to  determine  if  the  3-week  duration  of 
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aging/weathering  procedure  was  adequate.  The  diuation  of  this  study  was  18  weeks.  Nominal 
Mn  treatment  concentrations  included  0, 10, 18, 31, 54, 94, 164, 287,  and  503  mg  kg*’.  Samples 
from  each  treatment  concentration  were  analyzed  for  exchangeable  Mn  concentrations  at  3-week 
intervals  to  determine  if  increase  in  duration  of  aging/weadiering  procedure  beyond 
3  weeks  affects  exchangeable  Mn  concentrations  (directly  related  to  bioavailable  Mn). 

2.6  Chemical  Extraction  and  Anal^^es■ 

Soil  was  analyzed  for  total  metal  concentrations  following  USEPA  Method  200.8 
(USEPA,  1994)  using  inductively  coupled  plasma  mass  spectrometry  (ICP-MS).  Additional 
analysis  was  done  to  determine  exchangeable  Mn  fraction.  Exchangeable  Mn  was  extracted 
from  soil  using  0.05A/  CaCb  with  agitation  on  a  reciprocating  shaker  for  24  h.  All  reagents  used 
in  extraction  of  chemicals  fixim  soils  were  either  reagent  or  trace  metal  grade,  and  purified  water 
was  used  throughout  the  analytical  studies.  Glassware  was  washed  with  phosphate-fi%e 
detergent  followed  by  rinses  with  tap  water,  purified  water,  nitric  acid  1%  (v/v),  and  finally  again 
widi  purified  water.  Analyses  of  exchangeable  Mn  concentrations  were  conducted  using  a 
Peikin-Elmer  5100  PC  Atomic  Absorption  Spectrophotometer  equipped  with  an  AS-90 
autosampler. 

2.7  Toxicity  .^sessment. 

Two  earthworm  toxicity  tests  were  used  to  determine  eifiier  acute  or  chronic 
eni^ints.  The  acute  test  used  was  a  14-day  survival  test  adapted  fiom  Greene  et  at.  (1989). 

The  enc^int  of  this  test  is  number  of  adult  survivors.  The  chronic  test  used  in  this  study  was  a 
21-<by  reproduction  test  (ISO  1998;  adapted  fiom  Van  Gestel  et  al,  1989).  The  endpoint  of  this 
test  is  number  of  cocoons  produced.  Guidelines  for  these  assays  were  originally  developed  for 
use  with  artificial  soil  (USEPA  Standard  Artificial  Soil),  however  our  research  showed  that  these 
tests  could  also  be  successfully  conducted  using  natural  soils  (Kuperman  and  Simini,  2004). 

2.7.1  Eardiworm  Survival  Test 

2.7. 1 . 1  Principle  of  die  Test 

Adult  E.fetida  are  exposed  to  a  range  of  concentrations  of  the  test  chemical  added 
to  soil.  The  test  consists  of  two  steps:  first,  a  range-finding  test  in  which  adult  survival  is  assessed 
using  few  treatment  concentrations  (five)  and  replicates  (two);  and  second,  a  definitive  test  in  which 
survival,  Uve  wei^t,  and  dry  weight  are  assessed  using  a  greater  number  of  concentrations  and 
r^licates.  The  duration  of  each  test  is  2  weeks.  The  number  of  adult  survivors  in  treated  soils  is 
compared  to  the  number  in  die  controI(s)  to  quantify  ecotoxicological  parameters.  These 
parameters  include  the  No  Observed  Effect  Concentration  (NOEC),  the  Lowest  Observed  Effect 
Concentration  (LOEC)  and  the  effective  concentration  that  causes  an  x  percent  reduction  in  adults, 
i.e.,  ECx  (e.g.  EC20,  EC50). 
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2.7.1.2 


Validity  of  the  Test 


The  validity  criteria  are  included  in  the  test  as  part  of  the  Quality  Control 
procedures.  They  include  the  following  performance  parameters  for  the  negative  controls: 

(1)  The  mean  mortality  does  not  exceed  10%  at  the  end  of  the  range-finding  and 
definitive  tests; 

(2)  The  coefficient  pf  variation  for  the  mean  number  of  survivors  is  <30%  at  the 

end  of  the  test. 

2.7. 1 .3  Test  Conditions 

(1)  Earthworms  were  acclimated  for  48  h  in  die  test  soil.  Worms  with  fiilly 
developed  clitella  were  selected  for  uniformity  and  purged  overnight  on  moist  filter  paper. 

(2)  Worms  were  rinsed  twice  with  de-ionized  water,  blotted  on  paper  towels, 
weighed  on  an  analytical  balance,  and  placed  on  the  surface  of  soil  in  each  of  four  400-ml 
(9  cm  diam),  glass  containers  witii  screw  caps.  The  worms  were  selected  randomly  for 
placement  across  treatments. 

(3)  Plastic  film  was  stretched  over  the  top  of  the  containws  and  secured  with  the 
screw  caps.  Three  small  holes  were  made  in  the  wn^)  with  a  push-pin  to  allow  for  air  exchange. 

(4)  Worms  were  incubated  under  continuous  light  at  22±2°C  for  14  days. 

2.7. 1.4  En(4)oint  Determination 

At  the  end  of  the  study,  the  number  of  surviving  earthworms  in  each  beaker  were 
counted  and  recorded.  These  procedures  are  summarized  below: 

(1)  Soil  was  emptied  into  a  clean  dry  pan.  The  number  of  earthworm  survivors 
per  beaker  was  counted  and  recorded. 

(2)  Treatment  means  were  calculated  and  the  data  woe  analyzed  by  Analysis  of 
Variance  (ANOVA),  and  mean  separation  was  performed  using  Fischer’s  Least  Significant 
Difference  (LSD)  test  using  p<0.05  criterion. 

2.7.2  Earthworm  Reproduction  Test 

At  the  start  of  the  21 -day  chronic  reproduction  assay,  earthworms,  soil,  and 
beakers  were  prepared  as  described  above,  except  that  a  2  g  bolus  of  al&lfa  food  (see  Culturing 
Conditions)  was  added  to  each  beaker.  The  soil  was  then  hydrated  to  95%  of  SSL  water-holding 
edacity. 
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2.7.2. 1  Endpoint  Determination 

After  21  days,  surviving  earthworms  were  counted  and  weighed  as  described 
above.  Cocoons  were  recovered  by  gently  agitating  the  soil  on  a  l-mm  sieve  with  water  until 
only  the  cocoons  remained  on  the  surface  of  the  sieve.  The  number  of  cocoons  per  container 
(from  five  eartiiworms  per  treatment  replicate)  was  recorded. 

2.122  Validity  of  Test 

(1)  The  mean  mortality  does  not  exceed  10%  at  the  end  of  the  range-finding  and 
definitive  tests; 

(2)  The  coefficient  of  variation  for  the  mean  number  of  cocoons  is  <30%  at  the 

end  of  the  test. 

2.8  Selection  of  Metal  Concentrations  in  the  Soil. 

Preliminary  letiial  (survival)  and  sublethal  (cocoon  production)  toxicity  tests 
(range-finding  tests)  were  performed  with  each  metal  to  determine  the  range  of  concentrations  to 
use  in  tiie  definitive  toxicity  tests.  Five  concentrations  were  tised  for  each  of  the  metals  for  both 
survival  and  cocoon  production  assays.  The  range  of  concentrations  varied  according  to  the  type 
of  test  (adult  survival  or  reproduction). 

Data  fix>m  flie  range  finding  tests  were  used  to  select  the  chemical  compound 
containing  each  metal  that  was  s^ropriate  to  use  for  definitive  toxicity  tests  with  E.fetida,  and 
to  determine  treatment  concentrations  for  definitive  tests.  Additional  consideration  in  chemical 
form  selection  was  given  to  chemical  solubility  in  water  and  the  effect  that  each  chemical  form 
had  on  soil  pH.  Concentrations  selected  and  increment  factors  used  to  determine  treatment 
concentrations  for  definitive  tests  are  shown  in  Table  2. 


Table  2.  Nominal  Ba,  Be,  Mn,  and  Sb  concentrations  selected  for  definitive  toxicity  studies 
with  E.  fetida. 


Chemical 

Ba 

Be 

[n 

Sb 

Repro- 

Repro- 

Repro- 

Repro- 

Toxicity  Test 

Survival 

duction 

Survival 

duction 

Survival 

duction 

Survival 

duction 

Nomina] 

1  Concentration  (m 

g  kg'*) 

800 

112.5 

50 

20 

238 

100 

245 

60 

944 

225 

70 

28 

333.2 

274 

318 

86 

1114 

451 

98 

39 

478 

384 

412 

104 

1314 

519 

137 

55 

653 

537 

537 

124 

1551 

597 

192 

77 

914 

753 

617 

149 

1830 

686 

269 

108 

1278 

903 

697 

179 

2160 

789 

376 

151 

1792 

1054 

215 

2584 

907 

1043 
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Controls  included  positive  (30  mg  kg'‘  4-Nitroplienol),  negative  (no  chemicd 
added)  and  sulfate  (CaS04).  Sulfate  controls  were  based  on  estimated  sulfate  founts  in  hipest 
treatment  concentrations,  and  were  7,000  and  35,000  mg  kg'*  S04^'.  Four  replicates  were  used 
for  each  treatment  concentration  and  controls. 

2.9  Data  Analysis. 

Nonlinear  regression  procedures  were  applied  to  the  cocoon  production  data. 
Nonlinear  models  included  BCx  3s  a  parameter  to  determme  the  metal  concentration  producing  a 
specified  percentage  reduction  in  juvenile  production.  These  *‘x”  parameters  included  EC20 
EC50  levels.  The  first  parameter  (EC20)  is  the  benchmark  value  preferred  for  Eco-SSL 
development  for  soil  invertebrates,  and  the  second  (EC50),  more  commonly  used  in  the  past,  was 
included  to  enable  comparisons  of  the  results  produced  in  this  study  witii  results  reported  by 
other  researches.  The  asymptotic  standard  error  (a.s.e.)  and  95%  confidence  intervals  (C.I.) 
associated  with  the  point  estimates  woe  determined. 

Reproduction  data  were  analyzed  using  the  Gompertz  model  (Y- ax 
p)]x[c/ECp]Afr)^  described  in  Stephenson  et  al,  2001.  Histograms  of  the  residuals  and  stem-and-leaf 
graphs  were  examined  to  ensure  that  normality  assumptions  were  met,  aldiough  normality  is  not 
particularly  important  for  these  analyses  as  long  as  the  data  are  approximately  symmetrically 
distributed  (Stephenson  et  al,  2001).  Variances  of  the  residuals  were  examined  to  decide 
whether  or  not  to  weight  the  data,  and  to  select  potential  models.  The  Gompertz  model  had  the 
best  fit  for  all  metal  data.  The  fit  of  the  line  was  closest  to  the  data  points,  the  variances  were  the 
smallest  and  the  residuals  had  the  best  appearance  (i.e.,  most  random  scattering). 

Analysis  of  variance  (ANOVA)  was  used  to  determine  the  bounded  NOEC  and 
LOEC  values  in  both  range  finding  tests  and  definitive  tests.  These  analyses  woe  applied 
to  adult  survival  data,  and  to  cocoon  production  data.  Mean  separations  were  done  using 
Fisher’s  LSD  (Least  Significant  Difference)  pairwise  comparison  tests.  Significance  level 
oiP<  0.05  was  accepted  for  determining  the  NOEC  and  LOEC  values.  The  appropriate  tests  for 
assumptions  of  the  models  wrae  applied  as  described  above  before  using  either  ANOVA  or 
regression  procedures.  All  analyses  were  performed  using  SYSTAT  7.0.1  (SPSS,  Inc.,  1997). 

Raw  data  for  range-finding  and  definitive  tests  were  tabulated  and  are  listed  in 
Appendixes  A  and  B,  respectively.  Detailed  results  of  statistical  analysis  of  toxicity  test  data  are 
listed  in  Appradix  C. 


3.  RESULTS 

3.1  Soil  Analysis. 

Analysis  of  negative  control  soil  showed  that  Be  concentration  in  natural  SSL  soil 
used  in  this  study  was  below  the  mediod  the  detection  limit  (MDL)  of  2.5  mg  kg  *.  Total  Be 
concentrations  in  the  experimental  treatments  ranged  firom  95  to  124%  and  averaged  107%  of 
nominal  (Table  3). 
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The  natural  background  Mn  concentration  determined  in  the  negative  control 
treatment  was  94  mg  kg'*.  Total  extractable  Mn  concentrations  (in  excess 
of  background)  in  the  experimental  treatments  ranged  from  50  to  1 17%  and  averaged  94% 
of  nominal  (Table  3).  Exchangeable  Mn  fraction  expressed  as  percent  of  total  concentration 
increased  with  increasing  soil  Mn  loads  (Table  4).  There  were  no  trends  within  any  treatment 
concentration  in  the  amount  of  exchangeable  Mn  fraction  beyond  3  weeks  during  the  18-week 
aging/weathering  study.  These  results  confirmed  that  the  3-week  duration  for  simulated 
aging/weadiering  procedure  used  in  to  the  definitive  study  design  was  adequate  for  the  Eco-SSL 
benchmark  development. 

Analytical  procedures  for  Sb  determination  did  not  confirm  agreement  with  the 
nominal  treatment  concentrations.  Total  Sb  treatment  concentrations  determined  using  USEPA 
Method  200.8  ranged  fix)m  4  to  21%  and  averaged  8%  of  nominal  concentration.  These  results 
showed  that  this  standard  method  was  not  sufficient  for  total  Sb  analysis  in  SSL  soil.  Additional 
effort  was  made  to  improve  the  analytical  procedure.  Soils  were  digested  using  procedures 
described  in  SW-846  Method  3050B  (USEPA,  1996).  This  improved  the  efficiency  of  Sb 
extraction,  however  it  remained  relatively  low  and  averaged  58%  of  nominal  concentration.  For 
this  reason,  nominal  Sb  concentrations  were  used  in  determining  ecotoxicological  parameters  for 
Sb;  however  because  ERA  relies  on  the  determination  of  soil  concentrations  extracted  fiom  soil, 
toxicity  parameters  determined  finm  nominal  concentrations  may  have  to  be  adjusted  to  58%  of 
their  values  before  determining  an  Sb  Eco-SSL  to  best  conservatively-correspond  to  the  level  of 
Sb  extracted  fiom  soil  at  specific  levels  of  Sb  toxicity  in  soil. 

The  SSL  soil  pH  value  of  5.29  was  within  the  range  of  Eco-SSL’s  soil  matrix  of 
properties  tiiat  support  high  bioavailability  of  cationic  metals  in  natural  soils.  Soil  pH  decreased 
consistently  witii  increasing  chemical  loads  but  the  decrease  did  not  exceod  one  pH  unit  for  Ba, 
Mn,  and  Sb  treatments  (Table  5).  The  decrease  in  the  highest  Be  treatment  was  1.46  pH  units 
compared  with  untreated  SSL  soil  (negative  control).  In  the  sulfate  control,  soil  pH  decreased 
by  less  tiian  1 .0  pH  unit  in  both  7000  and  35000  mg  kg'*  S04^'  treatments  compared  with 
negative  control. 

3.2  Earthworm  Toxicity  Tests. 

3.2.1  Barium 

Range-finding  tests  for  barium  (Ba)  were  performed  using  Ba  nitrate  (BaN206). 
Nominal  concentrations  of  Ba  (w/w)  in  soil  were  0, 100, 500, 1000, 5000,  and  10000  mg  kg'*. 
Results  showed  cocoon  production  was  67%  of  control  at  100  mg  kg'*,  and  22%  of  control  at 
500  mg  kg'*.  No  cocoons  were  produced  above  500  mg  kg'*. 

ANOVA  results  of  definitive  cocoon  production  tests  using  actual  measured 
concentrations  of  Ba  (Table  3),  produced  a  bounded  NOEC  of 258  mg  kg'*  and  a  bounded  LOEC 
of  433  mg  kg'  at  P  <  0.05  (Table  6).  Cocoon  production  was  analyzed  by  nonlinear  regression. 
The  Gompertz  nonlinear  model  produced  best  fit  (Figure  1).  The  EC20  derived  fix)m  this  model 
for  toxicity  of  Ba  to  E.fetida  cocoon  production  was  370  mg  kg'*  and  the  EC50  was  664  mg  kg'* 
(Table  6). 
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Results  of  the  definitive  adult  survival  tests,  using  actual  measured  concentrations 
of  Ba  (Table  3),  showed  a  bounded  NOEC  of  1348  mg  kg*‘  and  a  bounded  LOEC  of  1585  mg  kg  ' 
(Table  6). 


Table  3.  Nominal  and  measured  concentrations  of  metals  in  soil. 


Chemical 

Ba  Concentration  (mg  kg"') 

Be  Concentra 

ion  (mg  kg’*) 

Survival 

Survival 

mssm 

nominal 

measured 

nominal 

measured 

nominal 

measured 

nominal 

measured 

inherent 

34 

inherent 

34 

inherent 

<2.5 

inherent 

<2.5 

800 

1000 

112.5 

153 

50 

79 

20 

24 

944 

1124 

225 

258 

70 

83 

28 

28 

1114 

1222 

451 

433 

98 

no 

39 

44 

1314 

1348 

519 

578 

137 

144 

55 

57 

1551 

1585 

597 

689 

192 

191 

77 

83 

1830 

2000 

686 

744 

269 

308 

108 

121 

2160 

2194 

789 

791 

376 

380 

151 

154 

2584 

2697 

907 

1000 

1043 

1222 

Mean  %  nominal 

109 

112 

116 

108 

Chemical 

Mn  Concentration  (mg  kg’‘) 

"Sb  Concentration  (mg  kg“*) 

Survival 

Survival 

■2S 

mssm 

nominal 

measured 

nominal 

measured 

nominal 

measured 

nominal 

measured 

inherent 

94 

inherent 

94 

inherent 

<2.5 

inherent 

<2.5 

238 

326 

100 

202 

245 

-- 

60 

— 

332 

449 

274 

386 

318 

— 

86 

— 

478 

611 

384 

528 

412 

— 

104 

653 

767 

537 

697 

537 

— 

124 

— 

914 

nil 

753 

1067 

617 

— 

149 

— ■ 

1278 

1444 

903 

nil 

697 

— 

179 

— 

1792 

2222 

1054 

1236 

215 

Mean  %  nominal 

125 

142 

’Several  extraction  methods  were  attempted  to  measure  Sb  levels  in  tiie  soil.  The  best  extraction  efficiency  was 
58%  of  the  nominal  Sb  concentration.  Nominal  Sb  concentrations  were  used  for  statistical  analyses. 
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Table  4.  Initial  total  manganese  concentrations  and  exchangeable  manganese  fractions  during 
18-week  aging/weathering  study  using  SSL  soil  amended  with  manganese  sulfate. 


Nominal 

Mn 

treatment 


Exchangeable  Mn  fraction  (%  of  total) 


Treatment 
mean 


Week  3 

Week  6 

Week  9 

Week  12 

Week  15 

Week  18 

(%  of  total) 

0 

5.4 

4.9 

7.3 

6.6 

6.2 

7.7 

6.4 

10 

18.0 

16.3 

19.9 

20.1 

16.3 

17.8 

18.1 

18 

27.1 

25.6 

28.7 

30.1 

23.5 

27.9 

27.2 

31 

42.3 

37.3 

39.1 

44.2 

38.8 

40.5 

40.4 

54 

60.1 

52.4 

54.9 

60.4 

48.5 

54.5 

55.1 

94 

85.8 

75.9 

76.0 

82.4 

65.3 

76.7 

77.0 

164 

75.2 

63.9 

66.7 

70.7 

56.3 

68.9 

66.9 

287 

106.3 

93.8 

94.3 

98.5 

82.2 

95.8 

95.2 

503 

127.3 

99.8 

104.7 

110.4 

101.7 

90.3 

105.7 

Table  5.  Summary  of  soil  pH  data.  Values  were  determined  following  a  3-week 
aging/weathering  procedure  in  studies  of  nominal  berylliiun,  manganese, 
antimony,  and  Ba  concentrations  amended  individrially  in  SSL  soil.  Data 
include  levels  used  in  both  survival  and  reproduction  tests  with  E.  fetida. 


Ba 

pH 

Be 

pH 

Mn 

pH 

Sb 

pH 

0 

5.29 

0 

5.29 

0 

5.29 

0 

5.29 

112.5 

5.00 

20 

4.89 

100 

5.11 

60 

5.14 

225 

4.88 

28 

4.75 

238 

4.92 

86 

5.09 

451 

4.72 

39 

4.66 

274 

4.88 

104 

5.08 

519 

4.66 

50 

4.55 

332 

4.85 

124 

5.01 

597 

4.63 

55 

4.67 

384 

4.87 

149 

4.96 

686 

4.63 

70 

4.36 

478 

4.81 

179 

4.91 

789 

4.54 

77 

4.37 

537 

4.81 

215 

4.82 

800 

4.54 

98 

4.26 

653 

4.76 

245 

4.82 

907 

4.50 

108 

4.26 

753 

4.74 

318 

4.69 

944 

4.50 

137 

4.13 

903 

4.72 

412 

4.66 

1043 

4.48 

151 

4.13 

914 

4.71 

537 

4.61 

1114 

4.45 

192 

4.02 

1054 

4.71 

617 

4.50 

1314 

4.44 

269 

[3.94 

1278 

4.68 

697 

4.39 

1551 

4.38 

376 

3.83 

1792 

4.58 

1830 

4.36 

2160 

4.29 

2584 

4.26 

3.2.2 


Beryllium 


Range-finding  tests  for  beryllium  (Be)  were  performed  using  Be  sulfate 
(BeS04*4H20).  Nominal  concentrations  of  Be  (w/w)  in  soil  were  0, 1, 10, 100, 500,  and 
1000  mg  kg'*.  Results  showed  100%  survival  rate  up  to  100  mg  kg'*.  All  adults  were  dead  at 
500  mg  kg°  and  1000  mg  kg'*. 

The  results  of  the  cocoon  tests  at  39  mg  kg'*  were  invalid  because  the  earthworms 
received  no  food  and  no  cocoons  were  produced.  These  data  were  therefore  not  included  in  the 
statistical  analyses.  ANOVA  results  of  definitive  cocoon  production  tests  using  actual  measured 
concentrations  of  Be  (Table  3)  produced  a  bounded  NOEC  of  57  mg  kg  *  and  a  bounded  LOEC 
of  83  mg  kg'*  at  P  <  0.05  (Table  6).  Cocoon  production  was  analyzed  by  nonlinear  regression. 
The  Gompertz  nonlinear  model  produced  best  fit  (Figure  1).  The  EC20  derived  fix)m  tiiis  model 
for  toxicity  of  Be  to  E,  fetida  cocoon  production  was  52  mg  kg'*  and  the  EC50  was  63  mg  kg' 
(Table  6). 


Results  of  the  definitive  adult  survival  tests  using  actual  measured  concraitrations 
of  Be  (Table  3)  showed  a  bounded  NOEC  of  83  mg  kg'*  and  a  bounded  LOEC  of  1 10  mg  kg'* 
at  P  <  0.05  (Table  6;  Figure  2). 

3.2.3  Manganese 

Range-finding  tests  for  manganese  (Mn)  were  performed  using  Mn  sulfate 
(MnS04»H201.  Nominal  concentrations  of  Mn  (w/w)  in  soil  were  0, 100, 500, 1000, 5000,  and 
10000  mg  kg'*.  Results  showed  cocoon  production  was  69%  of  control  at  500  mg  kg'  and  4% 
of  control  at  1000  mg  kg'*.  No  cocoons  were  produced  in  soil  containing  mg  more  Aan 
1000  mg  kg'*. 

ANOVA  results  of  definitive  cocoon  production  tests  using  actual  measured 
concentrations  of  Mn  ^able  3)  produced  a  bounded  NOEC  of  1 1 1 1  mg  kg  *  and  a  bounded 
LOEC  of  1236  mg  kg'*  at  P  <  0.05  (Table  6).  Cocoon  production  was  analyzed  by  nonlinear 
regression.  The  Gompertz  nonlinear  model  produced  best  fit  (Figure  1).  The  EC20  derived  from 
this  model  for  toxicity  of  Mn  to  E.  fetida  cocoon  production  was  629  mg  kg'*  (Table  6).  The 
EC50  was  927  mg  kg'^  (Table  6). 

Results  of  the  definitive  adult  survival  tests  using  actual  measured  concentrations 
of  Mn  (Table  3),  showed  a  bounded  NOEC  of  1444  mg  kg'*  and  a  bounded  LOEC  of 
2222  mg  kg'*  (Table  6).  Adult  EC20  and  EC50  values  were  1718  mg  kg'*  and  1920  mg  kg' , 
respectively  (Table  6). 
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Table  6.  Ecotoxicological  parameters  associated  with  exposure  of  E.  fetida  to  barium  (Ba), 
beryllium  (Be),  manganese  (Mn),  and  antimony  (Sb).* 


NOEC  =  no  observed  effects  concentration.  LOEC  =  lowest  observed  effect  concentration.  NOEC  and  LOEC 
were  derived  from  Analysis  of  Variance  and  comparison  of  mean  metal  concentrations  by  Fisher’s  LSD  test 
at  P  <  0.05.  EC20  =  effective  concentration  at  which  20  %  reduction  occurs.  EC50  "=  effective  concentration  at 
which  50  %  reduction  occurs.  Cocoon  EC20  and  EC50  were  derived  from  nonlinear  regression  using  the  Gompertz 
logistical  model  F  =  a  x  Ecotoxicological  parameters  for  Sb  are  based  on  nominal  Sb 

concentrations,  used  in  statistical  analyses. 
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Figure  1.  Concentration-response  curves  for  four  metals  vs.  cocoon  production  in  Eiseniafetida  reproduction 
tests  using  Gompertz  logistical  model  7  =  a  x  g(f'ogd-/’)]x[C/ECp]A6) 


Beryllium 


3.2.4 


Antimony 


Range-finding  tests  for  antimony  (Sb)  were  performed  using  Sb  sulfate  (SbS04). 
Nominal  concentrations  of  Sb  (w/w)  in  soil  were  0, 1, 10, 100, 500,  and  1000  mg  kg  ‘.  Results 
showed  100%  survival  rate  up  to  100  mg  kg'*.  All  adults  were  dead  at  500  mg  kg'*  and  1000  mg 
kg‘. 


ANOVA  results  of  definitive  cocoon  production  tests  using  actual  measured 
concentrations  of  Sb  (Table  3),  produced  a  bounded  NOEC  of  60  mg  kg'*  and  a  bounded  LOEC 
of  86  mg  kg'*  @  p  <  0.05  (Table  6).  Cocoon  production  was  analyzed  by  nonlinear  regression, 
based  on  nominal  concentration  values  of  Sb.  The  Gompertz  nonlinear  model  produced  best  fit 
(Figure  1).  The  EC20  derived  firom  this  model  for  toxicity  of  Sb  to  E.  fetida  cocoon  production 
was  30  mg  kg'*  and  the  EC50  was  70  mg  kg'*  (Table  6). 

Results  of  the  definitive  adult  survival  tests  using  actual  measured  concentrations 
of  Sb  (Table  3),  showed  a  bounded  NOEC  of  617  mg  kg'*  and  a  bounded  LOEC  of  697  mg  kg' 
(Table  6). 


4.  DISCUSSION 


Development  of  screening  level  benchmarks  for  Ecological  Risk  Assessment 
(ERA)  of  contaminated  soils  has  become  a  critical  need  in  recent  years  (USEPA,  2000).  To 
address  this  problem,  the  USEPA  in  conjunction  with  stakeholders  is  developing  Eco-SSLs  to 
identify  concentrations  of  chemicals  in  soil  that,  when  not  exceeded,  theoretically  protective 
of  terrestrial  ecosystems  within  specific  soil  bormdary  conditions  fi-om  unacceptable  harmful 
effects.  An  extensive  review  of  literature  (USEPA,  2000)  determined  that  there  was  insufficient 
information  for  Be,  Mn,  Sb,  and  Ba  to  generate  Eco-SSL  benchmarks  for  soil  invertebrates.  Our 
toxicity  studies  were  designed  to  specifically  fill  this  knowledge  gap. 

The  majority  of  soil  toxicity  tests  reported  in  the  literature  used  standard  artificial 
soil  with  high  organic  matter  content  (10%)  and  near  neutral  pH.  In  contrast,  we  select^  SSL 
soil  to  meet  the  criteria  for  Eco-SSL  development,  in  large  part  because  it  has  characteristics 
supporting  relatively  high  bioavailability  of  cationic  metals.  In  addition,  our  aging/weathering 
procedure  of  the  soils  loaded  with  the  range  of  metal  concentrations  allowed  us  to  more 
realistically  assess  the  toxicity  under  conditions  more  closely  resembling  the  potential  toxic 
effects  of  Be,  Mn,  Sb,  and  Ba  in  the  field.  This  study  was  designed  to  produce  benchmark  data 
for  use  in  developing  Eco-SSLs  for  soil  invertebrates  for  Ba,  Be,  Mn,  and  Sb.  Thus  results  firom 
our  study  may  not  directly  compare  to  those  of  ofiier  studies  in  the  literature,  since  none  of  them 
were  designed  to  specifically  quantify  metal  toxicity  to  soil  invertebrates  using  the  Eco-SSL 
requirements  for  toxicity  testing  of  metals. 


The  present  study  has  determined  soil  toxicity  threshold  parameters  NOEC, 
LOEC,  EC20,  and  EC50  for  earthworms  in  response  to  soil  contamination  by  Ba,  Be,  Mn,  or  Sb, 
respectively.  Earthworm  is  considered  to  be  an  important  component  of  the  ecosystem  for 
metabolizing  organic  matter,  aerating  the  soil,  and  serving  as  food  for  higher  organisms  in  food 
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webs.  Several  species  of  earthworm  have  been  used  by  investigators  to  assess  toxicity  caused  by 
a  number  of  metal  contaminants  in  soils.  Ecotoxicological  parameters  for  earthworms  have  been 
well  characterized  for  cadmium,  copper,  lead,  and  zinc  (Kula  and  Larink,  1998;  Spurgeon,  et  al, 
1994,  Van  Gestel,  et  al.,  1989).  However,  little  work  has  been  done  to  characterize  earthworm 
soil  toxicity  thresholds  for  Ba,  Be,  Mn,  and  Sb  in  either  artificial  or  natural  soils. 

Our  results  indicate  that  Sb  and  Be  are  the  most  toxic  to  E.  fetida  of  the  four 
metals  tested  in  this  study.  Cocoon  production  EC20  for  Sb  (nominal)  and  Be  were  30  mg  kg'* 
and  52  mg  kg  ,  respectively  (Table  6).  Comparatively,  toxicity  of  Ba  and  Mn  to  E.  fetida  was 
much  less.  Cocoon  production  EC20  for  Ba  and  Mn  were  370  mg  kg'*  and  629  mg  kg'*, 
respectively  (Table  6).  Reproductive  endpoints  in  all  tests  were  more  sensitive  compared  with 
adult  survival  (Tables  6).  Spurgeon  et  al,  (1994)  found  a  similar  response  of  E.  fetida  to 
cadmium,  copper,  lead,  and  zinc.  Cocoon  production  EC50  values  were  46.3  mg  kg'*, 

53.3  mg  kg  *,  1,940  mg  kg  *,  and  276  mg  kg'*  for  Cd,  Cu,  Pb,  and  Zn,  respectively.  Adult 
survival  EC50  in  the  Spurgeon  et  al.  (1994)  study  ranged  from  2  to  13  times  greater  than  cocoon 
production  EC50.  We  found  similar  results.  Adult  survival  LOEC  in  our  study  ranged  from  2.5 
to  8  times  greater  than  cocoon  production  LOEC.  These  results  support  the  Eco-SSL 
requirement  of  the  use  of  reproductive  endpoints  for  benchmark  development. 

Beryllium  is  one  of  tiie  least  studied  metals  regarding  its  effects  on  soil 
invertebrates,  although  it  is  considered  one  of  the  problem  metals  of  the  future  (Newland,  1982). 
It  is  a  component  of  various  fossil  fuel  types  and  is  increasingly  used  in  aircraft  industry,  space 
research,  nuclear  energy  development  (Ireland,  1986),  X-ray  tube  windows  manufacturing  and  in 
production  of  non-sparking  tools  composed  of  copper-beryllium  alloy  (Thorat  et  al,  2001). 
Beryllium  concentrations  in  Aberdeen  Proving  Cirround  soil  (including  contaminated  sites)  in  the 
areas  adjacent  to  soil  collection  ranged  from  0.3  to  1.4  mg  kg'*  (Hlohowskyj  et  al,  1999). 
Extensive  toxicological  studies  of  exposure  effects  in  humans  and  experimental  animals  have 
established  that  Be  can  cause  pulmonary  and  systemic  granulomatous  disease  known  as  chronic 
Be  ^sease  (Sprince  and  Kazami,  1980),  necrosis  and  tumors  in  animals  (Witschi,  1971),  can 
mhibit  certain  enzymes,  including  alkaline  phosphatase  (Reiner,  1971),  and  can  inhibit  plant  and 
animal  growth  (Newland,  1982).  Ireland  (1986)  reported  increased  mortality  and  growth 
suppression  in  a  terrestrial  snail  Achatina  fiilica  (Pulmonata)  fed  10  pg  ml'*  Be  in  the  diet 
containing  the  sub-optimal  calcium  concentrations.  Beryllium,  along  with  Sb,  was  the  most 
toxic  metal  among  the  four  chemicals  tested  in  our  study,  and  the  estimated  ecotoxicological 
parameters  for  E.  fetida  are  the  first  in  the  available  literature  for  a  soil  invertebrate  species. 

Few  studies  have  investigated  Sb  concentrations  in  soil  (Cal -Prieto  et  al,  2001; 
Crecelius  et  al,  1974;  Kabata-Pendias  and  Pendias,  1992;  van  der  Voet  and  de  Wolff,  1996). 
Reported  concentrations  ranged  from  0.17  mg  kg'*  in  organic  soils  in  Norway  to  1489  mg  kg'*  in 
vicinity  of  Sb  smelter  in  northeast  England  (Ainsworth  and  Cooke,  1991),  and  corresponded  with 
treatment  concentrations  used  in  our  study.  Antimony  concentrations  in  soil  (including 
contaminated  sites)  at  the  Aberdeen  Proving  Ground  in  the  areas  adjacent  to  the 
location  where  the  SSL  soil  was  collected  ranged  from  0.1  to  501  mg  kg'*  (Hlohowskyj  et  al, 
1999).  No  information  could  be  found  in  the  available  literature  on  ecotoxicological  effects  of 
Sb  to  soil  invertebrates.  Developing  such  information  is  especially  important  since  input  to  the 
soil  ecosystems  was  estimated  at  26000 1  y*  of  Sb  (Cal-Prieto  et  al,  2001).  This  anthropogenic 
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contribution  of  Sb  is  10-fold  higher  compared  with  the  Sb  emissions  from  natural  sources  (ca. 
2600 1  y  *)  reported  by  Nriagu  (1990).  Limited  data  for  soil  biota  were  reported  by  Rafel  and 
Popov  (1 988)  as  part  of  validation  efforts  for  developing  the  USSR  maximum  allowable 
concentrations  of  Sb  in  soil.  These  authors  reported  23-52%  reduction  in  seed  germination  and 
26-62%  reduction  in  root  growth  at  1002  mg  kg'‘  Sb  in  tests  with  barley,  wheat,  radish,  pees,  and 
onion.  Decrease  in  ammonia  mineralization  and  nitrate  accumulation  was  observed  at  Sb 
concentrations  of  52  and  102  mg  kg'*  in  their  study.  Other  measures  of  soil  biological  activity 
were  also  affected,  including  decrease  in  soil  enzyme  catalase  activity  and  stimulation  of  soil 
respiration  at  102  mg  Sb  kg'*  (Rafel  and  Popov,  1988). 

Difficulties  encoimtered  with  the  efficiency  of  extraction  of  Sb  that  is 
aged/weathered  in  soil  prior  to  analytical  determination,  using  natural  SSL  amended  with  Sb, 
may  be  symptomatic  of  a  larger  problem  regarding  chemical  characterization  data  during  ERA 
activities  at  contaminated  sites.  Low  Sb  recovery  rates  using  standard  USEPA  methods  suggest 
that  trae  concentrations  of  this  metal  will  be  underestimated  during  site  characterization  efforts. 
The  recovery  rates  of  8  and  58  %  determined  for  Sb  aged/weathered  in  soil  in  our  study,  using 
USEPA  mefliods  200.8  and  3050B  respectively,  were  below  recovery  rates  of  70  and  88  % 
preAdously  reported  for  freshly-spiked  soils.  This  clearly  indicates  that  USEPA  method  3050B 
appears  better  suited  to  extract  aged/weathered  Sb  from  soil ,  such  as  that  which  typically  occurs 
at  Superfund  and  other  contaminated  sites,  and  this  potential  discrepancy  in  extractability  should 
be  corrected  for  at  the  time  of  compilation  of  a  list  of  contaminants  of  potential  ecological 
concern  (COPEC)  in  the  screening  phase  of  ERA.  To  use  the  ecotoxicological  parameters  from 
this  study,  which  are  based  on  nominal  Sb  values,  it  is  recommended  that  these  nominal  Sb 
values  be  adjusted  to  58%  of  nominal  to  account  for  the  aging/weathering  of  Sb  in  soil  (i.e., 
adjusted  to  58%  of  nominal  prior  to  determining  the  Eco-SSL).  Aging/weathering  of  Sb  in  soils 
typically  occurs  even  more  extensively  in  the  field,  but  simulated  aging/weathering  provides  a 
conservative  estimate  of  what  might  otherwise  be  extractable  from  field  soils.  This  is  especially 
important  given  the  steep  slope  of  the  concentration-response  curve  for  reproductive  endpoint 
determined  from  the  Earthworm  Reproduction  Test  in  our  study  (Figure  1.),  which  establishes  a 
narrow  toxicity  threshold  range  from  30  to  70  mg  kg'*  based  on  EC20  and  EC50  estimates, 
respectively  (Table  6).  The  43  %  difference  between  these  two  estimates  is  within  the  potential 
recovery  error  rate  of  analytical  methods  used.  Disregarding  this  potential  error,  especially 
without  adjustment  of  the  Eco-SSL  for  aging/weathering,  can  ofiierwise  lead  to  a  removal  of  Sb 
from  the  COPEC  list  while  its  extracted  concentrations  represent  field  concentrations  toxic  to 
relevant  ecological  receptors.  Adjustment  of  the  values  of  the  ecotoxicological  parameters 
determined  from  nominal  concentrations,  prior  to  determination  of  the  Eco-SSL,  is  properly  left 
to  those  evaluating  benchmarks  for  Eco-SSL  development;  however,  in  these  studies  an 
adjustment  to  58%  of  nominal  corresponds  to  the  mean  recovery  rate  following  3  weeks  of 
aging/weathering  of  Sb  in  soil. 

Natural  Mn  concentration  in  SSL  soil  of  94  mg  kg'*  was  within  the  range  of  Mn 
concentrations  reported  for  soils  (including  contaminated  sites)  at  the  Aberdeen  Proving  Ground, 
which  ranged  fix)m  4.9  to  1 140  mg  kg'*  (Hlohowskyj  et  al,  1999).  Manganese  is  a  required 
nutrient  essential  for  plants  and  animals.  Manganese  was  the  most  previously  investigated  of  the 
four  metals  in  diis  study,  however  none  of  the  previous  studies  involved  invertebrate  exposures 
in  natural  soils.  Reinecke  and  Reinecke  (1996)  reported  reduction  in  growth  and  development 
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(measured  as  time  needed  for  clitellum  development)  of  E.fetida  fed  with  cattle  manure  spiked 
with  Mn  at  151.7  mg  kg  *,  This  value  falls  well  below  the  EC20  of  629  mg  kg'*  determined  in  our 
study.  In  a  later  study,  Reinecke  and  Reinecke  (1997)  reported  damage  to  spermatozoan 
structure  from  treatments  containing  food  spiked  with  Mn  at  61.57  mg  kg'*.  Nottrot  et  al  (1987) 
reported  no  effect  on  feeding  activity  and  growth  of  collembolan  Orchesella  cincta  fed  with 
green  algae  spiked  with  up  to  25  :mol  Mn  g  *  dry  mass,  however  that  study  was  conducted  on 
dental  plaster.  Joosse  and  van  Vliet  (1984)  reported  no  effect  on  respiration  of  woodlice  fed  with 
litter  containing  Mn  at  1000  mg  kg  *  on  a  porous  tile.  There  was  no  soil  exposure  incorporated 
in  that  study. 


Natural  Ba  concentration  in  SSL  soil  of  34  mg  kg'*  was  within  the  Ba 
concentrations  found  in  soils  (including  contaminated  sites)  at  the  Aberdeen  Proving  Ground, 
which  ranged  from  9.8  to  1580  mg  kg'*  (Hlohowsk>j  et  al,  1999).  Limited  Ba  ecotoxicological 
information  for  soil  invertebrates  is  available  from  literature.  Grace  (1990)  investigated  oral 
toxicity  of  Ba  metaborate  to  the  Eastern  Subterranean  Termite  Reticulitermes  flavipes  (Kollar)  in 
no-choice  assays  by  feeding  termite  workers  for  15  days  on  filter  papers  treated  with 
concentrations  of  Ba.  Results  of  this  study  are  similar  to  results  of  our  14-day  adult  survival 
defimtive  test.  Grace  (1990)  reported  19%  mortality  in  1780  mg  Ba  kg'*  treatment,  which  was 
comparable  with  24%  adult  mortality  at  1551  mg  Ba  kg'*  treatment  observed  in  our 
investigation.  However,  direct  comparisons  of  feeding  assay  results  with  soil  exposure  studies 
using  different  species  should  be  treated  with  caution. 


5.  CONCLUSIONS 

This  study  has  produced  ecotoxicological  benchmark  data  for  barium  (Ba), 
beryllium  (Be),  manganese  (Mn),  and  antimony  (Sb),  based  on  soil  toxicity  to  the  earthworm 
E.fetida,  for  use  in  the  developing  Eco-SSLs.  The  relative  toxicity  of  the  4  metals  tested  in  this 
study  was  Sb  (nominal)  >  Be  >  Ba  >  Mn.  When  the  ecotoxicological  parameters  for  Sb  are 
adjusted  by  58%  to  account  for  reduced  extractability  of  Sb  after  3  weeks  of  aging/weathering  in 
soil,  the  relative  toxicity  becomes  Sb  >Be  >  Ba  >  Mn.  It  is  strongly  recommended  that  the 
nominal  Sb  benchmark  values  from  this  study  be  adjusted  to  58%  of  nominal,  to  account  for  the 
aging/weathering  of  Sb  in  soil  (i.e.,  adjusted  to  58%  of  nominal  prior  to  determining  the  Eco- 
SSL).  Be  and  Sb  (nominal)  were  eight  to  ten  times  more  toxic  to  cocoon  production  in 
E.fetida  than  were  Ba  and  Mn.  Adult  survival  was  2.5  to  8  times  greater  than  cocoon  production 
indicating  that  reproduction  tests  provide  a  more  sensitive  evaluation  of  effect  than  survival  and 
therefore  should  be  used  to  set  screening  criteria.  These  tests  were  performed  using  a  natural 
soil,  Sassafias  sandy  loam.  Sassafias  sandy  loam  has  relatively  low  pH,  low  organic  matter,  low 
cation  exchange  capacity,  and  high  sand  content.  Such  characteristics  support  relatively  high 
bioavailability  of  cationic  metals  in  soil.  Furthermore,  aging  and  weathering  of  die  soil  produced 
a  soil  microenvironment  more  similar  to  field  conditions  than  previous  studies  where  soil 
invertebrates  were  exposed  immediately  following  spiking  of  soil.  These  study  results  will  be 
provided  to  the  Ecological  Soil  Screening  Level  (Eco-SSL)  workgroup  for  review.  Results  will 
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undergo  qiiality  control  review  by  die  Eco-SSL  task  group  before  inclusion  in  the  Eco-SSL 
database,  and  before  being  used  for  developing  Ecological  Soil  Screening  Levels  (Eco-SSLs)  for 
Be,  Mn,  Sb,  and  Ba. 
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APPENDIX  A 


RANGE-FINDING  TEST  DATA 


Range-finding  invertebrate  assays 
Fresh  SSL  soil 
Compound:  Ba  [BaS04] 

Start  Date:  19-Apr-OO 

Invertebrate:  E.  fetida 

Treatment  Rep  Adults  MEAN  Cocoons  MEAN  Reduction 


>pm 

4/10/00 

S.E. 

4/10/00 

S.E. 

% 

0 

1 

5 

0.0 

5 

5 

0 

0 

2 

5 

0.0 

5 

0.0 

100 

3 

5 

0.0 

1 

2.0 

60 

100 

4 

5 

0.0 

3 

1.0 

500 

1 

5 

0.0 

3 

4.5 

10 

500 

2 

5 

0.0 

6 

1.5 

1000 

1 

5 

0.0 

1 

1.5 

70 

1000 

2 

5 

0.0 

2 

0.5 

5000 

1 

5 

0.0 

3 

4.5 

10 

5000 

2 

5 

0.0 

6 

1.5 

10000 

1 

5 

0.0 

2 

2.0 

60 

10000 

2 

5 

0.0 

2 

0.0 

Range-finding  invertebrate  assays 
Fresh  SSL  soil 
Compound:  Ba  [BaO] 

Start  Date:  26-Jun-OO 

Invertebrate:  E.  fetida 


Treatment  Rep 

Adults 

MEAN  Cocoons 

MEAN 

Reduction 

Be  ppm 

7/17/00 

S.E. 

7/17/00 

S.E. 

% 

0 

1 

5 

5.0 

2.0 

2.0 

0 

0 

2 

5 

0.0 

2.0 

0.0 

100 

3 

5 

4.0 

1 

1.0 

50 

100 

4 

3 

1.0 

1 

0.0 

500 

1 

0 

0.0 

0 

0.0 

100 

500 

2 

0 

0.0 

0 

0.0 

1000 

1 

0 

0.0 

0 

0.0 

100 

1000 

2 

0 

0.0 

0 

0.0 
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Range-finding  invertebrate  assays 
Fresh  SSL  soil 

Compound:  Ba  [Ba(N03)2] 

Start  Date:  27 -Sep— 00 

Invertebrate:  E.  fetida 


Treatment  Rep 

Adults 

MEAN 

Cocoons 

MEAN 

Reduction 

Be  ppm 

10/18/00 

S.E. 

7/17/00 

S.E. 

% 

0 

1 

5 

5.0 

7 

4.5 

0 

2 

5 

0.0 

2 

2.5 

100 

1 

5 

5.0 

1 

3.0 

33 

100 

2 

5 

0.0 

5 

2.0 

500 

1 

5 

5.0 

0 

1.0 

78 

500 

2 

5 

0.0 

2 

1.0 

1000 

1 

5 

5.0 

0 

0.0 

100 

1000 

2 

5 

0.0 

0 

0.0 

5000 

1 

0 

0.0 

0 

0.0 

100 

5000 

2 

0 

0.0 

0 

0.0 

10000 

1 

0 

0.0 

0 

0.0 

100 

10000 

2 

0 

0.0 

0 

0.0 

Range-finding  invertebrate  assays 
Fresh  SSL  soil 


Compound:  Ba  [Ba(C2H302)2] 

Start  Date:  27-Sep-OO 

Invertebrate:  E.  fetida 


Treatment  Rep 

Adults 

MEAN 

Cocoons 

MEAN 

Reduction 

Be  ppm 

10/18/00 

S.E. 

7/17/00 

S.E. 

% 

0 

1 

4 

4.5 

3 

3.0 

0 

2 

5 

0.5 

3 

0.0 

100 

1 

5 

5.0 

1 

2.5 

17 

100 

2 

5 

0.0 

4 

1.5 

500 

1 

5 

5.0 

0 

0.0 

100 

500 

2 

5 

0.0 

0 

0.0 

1000 

1 

0 

0.0 

0 

0.0 

100 

1000 

2 

0 

0.0 

0 

0.0 

5000 

1 

0 

0.0 

0 

0.0 

100 

5000 

2 

0 

0.0 

0 

0.0 

10000 

1 

0 

0.0 

0 

0.0 

100 

10000 

2 

0 

0.0 

0 

0.0 
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Range-finding  invertebrate  assays 
Fresh  SSL  soil 
Compound:  Be  [BeS04] 

Start  Date :  1 7 -Apr-00 

Invertebrate:  E.  fetida 


Treatment  Rep 

Adults 

MEAN 

Cocoons 

MEAN 

Reduction 

Be  ppm 

5/10/00 

S.E. 

5/10/00 

S.E. 

% 

0 

1 

5 

5.0 

7 

6.0 

0 

2 

5 

0.0 

5 

1.0 

1 

1 

5 

5.0 

5 

4.5 

25 

1 

2 

5 

0.0 

4 

1.5 

10 

1 

5 

5.0 

8 

6.0 

0 

10 

2 

4 

4.5 

4 

2.0 

100 

1 

0 

1.0 

0 

0.0 

100 

100 

2 

2 

1.0 

0 

0.0 

500 

1 

0 

0.0 

0 

0.0 

100 

500 

2 

0 

0.0 

0 

0.0 

1000 

1 

0 

0.0 

0 

0.0 

100 

1000 

2 

0 

0.0 

0 

0.0 

Range-finding  invertebrate  assays 
Fresh  SSL  soil 
Compound:  Mn  [MnSO^ 

Start  Date:  19-Apr-OO 

Invertebrate:  E.  fetida 


Treatment  Rep 

Adults 

MEAN 

Cocoons 

MEAN 

Reduction 

Be  ppm 

5/10/00 

S.E. 

5/10/00 

S.E. 

% 

0 

1 

5 

4.0 

24 

13.0 

0 

2 

4 

0.0 

2 

10.0 

100 

1 

5 

4.0 

17 

14.0 

0 

100 

2 

4 

0.0 

11 

2.5 

500 

1 

5 

5.0 

14 

9.0 

31 

500 

2 

4 

4.5 

4 

5.0 

1000 

1 

0 

1.0 

0 

0.0 

96 

1000 

2 

2 

1.0 

1 

0.0 

5000 

1 

0 

0.0 

0 

0.0 

100 

5000 

2 

0 

0.0 

0 

0.0 

10000 

1 

0 

0.0 

0 

0.0 

100 

10000 

2 

0 

0.0 

0 

0.0 
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Range-finding  invertebrate  assays 
Fresh  SSL  soil 

Compound:  Sb  [Sb-d-tartrate] 

Start  Date:  17-Apr-OO 

invertebrate:  E.  fetida 


Treatment  Rep 

Adults 

MEAN 

Cocoons 

MEAN 

Reduction 

Be  ppm 

5/10/00 

S.E. 

5/10/00 

S.E. 

% 

0 

1 

5 

4.5 

4 

3.5 

0 

2 

4 

0.0 

3 

0.5 

1 

1 

5 

4.5 

6 

4.0 

0 

1 

2 

4 

0.0 

2 

2.5 

10 

1 

5 

5.0 

11 

7.5 

0 

10 

2 

4 

4.5 

4 

3.5 

100 

1 

4 

4.5 

3 

1.5 

58 

100 

2 

5 

0.5 

0 

1.5 

500 

1 

0 

0.0 

0 

0.0 

100 

500 

2 

0 

0.0 

0 

0.0 

1000 

1 

0 

0.0 

0 

0.0 

100 

1000 

2 

0 

0.0 

0 

0.0 
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APPENDIX  B 


DEFINITIVE  TEST  DATA 


Definitive  invertebrate  assays 
Aged  SSL  soil 

Compound:  Ba  [Ba(N03)2] 

Start  Date:  28-Nov-OO 

Invertebrate:  £.  fetida 
Survival 


Ba  (ppm) 

0 

0 

0 

0 

800 

800 

800 

800 

944 

944 

944 

944 

1114 

1114 

1114 

1114 

1314 

1314 

1314 

1314 

1551 

1551 

1551 

1551 

1830 

1830 

1830 

1830 

2160 

2160 

2160 

2160 

2548 

2548 

2548 

2548 


Rep  12/14 
Adults 

1 

2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 


Mean  S.E. 

5 

5 

5 

5  5.0  0.0 

5 
5 
5 

5  5.0  0.0 

5 
5 
5 

5  5.0  0.0 

5 
5 
5 

5  5.0  0.0 

5 
5 
5 

5  5.0  0.0 

5 
5 
5 

0  3.75  1.3 

0 
0 
0 

0  0.0  0.0 

0 
0 
0 

0  0.0  0.0 

0 
0 
0 
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Definitive  invertebrate  assays 
Aged  SSL  soil 

Compound:  Ba  [Ba(N03)2] 

Start  Date:  21  'Nov-OO 

Invertebrate:  E.  fetida  Reproduction 

Treatment  Rep  12/14 

Ba  (ppm)  Cocoons  Mean  S.E. 

0  1  27  21.0  2.0 

0  2  18 

0  3  19 

0  4  20 

112.5  1  22  21.8  0.6 

112.5  2  20 

112.5  3  22 

112.5  4  23 

225  1  20  17.0  2.0 

225  2  11 

225  3  18 

225  4  19 

451  1  16  15.0  0.9 

451  2  17 

451  3  13 

451  4  14 

519  1  11  14.0  2.3 

519  2  15 

519  3  10 

519  4  20 

597  1  17  12.3  2.4 

597  2  15 

597  3  6 

597  4  11 

686  1  9  8.5  0.3 

686  2  8 

686  3  9 

686  4  8 

789  1  5  5.5  1.3 

789  2  3 

789  3  5 

789  4  9 

907  1  2  5.0  1.2 

907  2  5 

907  3  5 

907  4  8 

1043  1  1  1.5  0.5 

1043  2  1 

1043  3  3 

1043  4  1 
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Definitive  invertebrate  assays 
Aged  SSL  soil 

Compound:  Be  [BeS04] 

Start  Date:  24-Aug-OO 

Invertebrate:  E.  fetida 
Survival 

Treatment  Rep  09/06 


Be  (ppm) 

Adults 

Mean 

S.E. 

0 

1 

5 

5.0 

0.0 

0 

2 

5 

0 

3 

5 

0 

4 

5 

50 

1 

5 

5.0 

0.0 

50 

2 

5 

50 

3 

5 

50 

4 

5 

70 

1 

4 

4.8 

0.3 

70 

2 

5 

70 

3 

5 

70 

4 

5 

98 

1 

5 

3.75 

0.5 

98 

2 

3 

98 

3 

3 

98 

4 

4 

137 

1 

5 

4.5 

0.3 

137 

2 

5 

137 

3 

4 

137 

4 

4 

192 

1 

0 

1.0 

0.4 

192 

2 

2 

192 

3 

1 

192 

4 

1 

269 

1 

0 

0.0 

0.0 

269 

2 

0 

269 

3 

0 

269 

4 

0 

376 

1 

0 

0.0 

0.0 

376 

2 

0 

376 

3 

0 

376 

4 

0 

PNP-30 

1 

0 

PNP-30 

2 

0 

PNP-30 

3 

0 

PNP-30 

4 

0 
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Definitive  invertebrate  assays 
Aged  SSL  soil 
Compound:  Be  [BeS04] 

Start  Date:  22*Aug-00 

Invertebrate:  E.  fetida 
Reproduction 

Treatment  Rep  09/13 


Be  (ppm) 

Cocoons 

Mean 

S.E. 

0 

1 

8 

7.3 

1.7 

0 

2 

11 

0 

3 

7 

0 

4 

3 

20 

1 

6 

8.0 

1.7 

20 

2 

7 

20 

3 

6 

20 

4 

13 

28 

1 

6 

6.8 

1.1 

28 

2 

10 

28 

3 

6 

28 

4 

5 

39 

1 

0 

0.0 

0.0 

39 

2 

0 

39 

3 

0 

39 

4 

0 

55 

1 

5 

5.0 

1.2 

55 

2 

5 

55 

3 

8 

55 

4 

2 

77 

1 

0 

0.3 

0.3 

77 

2 

0 

77 

3 

0 

77 

4 

1 

108 

1 

1 

1.0 

0.7 

108 

2 

0 

108 

3 

3 

108 

4 

0 

151 

1 

0 

0.0 

0.0 

151 

2 

0 

151 

3 

0 

151 

4 

0 

PNP-30 

1 

0 

PNP-30 

2 

0 

PNP-30 

3 

0 

PNP-30 

4 

0 
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Definitive  invertebrate  assays 
Aged  SSL  soil 

Compound:  Mn  [MnS04] 

Start  Date:  24-Aug-OO 

Invertebrate:  E.  fetida 
Survival 

Treatment  Rep  09/06 
Mn  (ppm)  Adults  Mean  S.E. 

0  1  5  5.0  0.0 

0  2  5 

0  3  5 

0  4  5 

238  1  5  5.0  0.0 

238  2  5 

238  3  5 

238  4  5 

332  1  5  5.0  0.0 

332  2  5 

332  3  5 

332  4  5 

478  1  5  5.0  0.0 

478  2  5 

478  3  5 

478  4  5 

653  1  5  5.0  0.0 

653  2  5 

653  3  5 

653  4  5 

914  1  5  4.8  0.5 

914  2  4 

914  3  5 

914  4  5 

1278  1  4  4.8  0.5 

1278  2  5 

1278  3  5 

1278  4  5 

1792  1  1  0.8  0.5 

1792  2  1 

1792  3  1 

1792  4  0 

1  0 

2  0 

3  0 

4  0 


PNP-30 

PNP-30 

PNP-30 

PNP-30 
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Definitive  invertebrate  assays 
Aged  SSL  soil 

Compound:  Mn  [MnS04] 

Start  Date:  22-Aug-OO 

Invertebrate:  E.  fetida 
Reproduction 

Treatment  Rep  09/13 
Mn(ppm)  Cocoons  Mean  S.E. 


0  1 

0  2 

0  3 

0  4 

100  1 

100  2 

100  3 

100  4 

274  1 

274  2 

274  3 

274  4 

384  1 

384  2 

384  3 

384  4 

537  1 

537  2 

537  3 

537  4 

753  1 

753  2 

753  3 

753  4 

903  1 

903  2 

903  3 

903  4 

1054  1 

1054  2 

1054  3 


1054  4 

PNP-30  1 

PNP-30  2 

PNP-30  3 

PNP-30  4 


11  7.0  3.7 

5 
3 

9 

17  13.8  2.9 

14 

10 
14 

14  8.8  5.1 

12 
3 

6 

9  7.3  4.3 

6 
12 
2 

6  6.3  4.7 

13 
3 
3 

3  3.3  1.7 

5 
1 

4 

2  2.8  2.2 

6 
2 
1 

3  1.5  1.0 

1 
1 
1 
0 
0 
0 
0 
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Definitive  invertebrate  assays 
Aged  SSL  soil 

Compound:  Sb  [Sb2(S04)3]] 

Start  Date:  24-Aug-OO 

Invertebrate:  E.  fetida 
Survival 

Treatment  Rep  09/06 
Sb  (ppm)  Adults  Mean  S.E. 

0  1  5  5.0  0.0 

0  2  5 

0  3  5 

0  4  5 

0  1  5  5.0  0.0 

0  2  5 

0  3  5 

0  4  5 

245  1  5  5.0  0.0 

245  2  5 

245  3  5 

245  4  5 

318  1  5  4.8  0.5 

318  2  5 

318  3  5 

318  4  5 

412  1  4  4.8  0.5 

412  2  5 

412  3  5 

412  4  5 

537  1  5  4.0  1.2 

537  2  4 

537  3  5 

537  4  5 

617  1  3  4.0  1.2 

617  2  5 

617  3  5 

617  4  3 

697  -1  4  3.5  1.3 

697  2  5 

697  3  3 

697  4  2 


1  0 

2  0 

3  0 

4  0 


PNP-30 

PNP-30 

PNP-30 

PNP-30 
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Definitive  invertebrate  assays 
Aged  SSL  soil 

Sb 

Compound:  [Sb2(S04)3]]] 

Start  Date:  22-Aug-OO 

Invertebrate:  E.  fetida 
Reproduction 

Treatment  Rep  09/11 

Sb  (ppm)  Cocoons  Mean  S.E. 

0  1  11  7.1  3.2 

0  2  5 

0  3  3 

0  4  9 

0  1  8 

0  2  11 

0  3  7 

0  4  3 

60  1  5  5.3  0.6 

60  2  6 

60  3  5 

60  4 

86  1  4  2.5  1.7 

86  2  4 

86  3  1 

86  4  1 

104  1  0  2.8  2.2 

104  2  5 

104  3  4 

104  4  2 

124  1  3  3.0  0.8 

124  2  2 

124  3  3 

124  4  4 

149  1  2  1.3  1.0 

149  2  0 

149  3  1 

149  4  2 

179  1  0  0.3  0.5 

179  2  0 

179  3  0 

179  4  1 

215  1  0  0.0  0.0 

215  2  0 

215  3  0 

215  4  0 

PNP-30  1  0 

PNP-30  2  0 

PNP-30  3  0 

PNP-30  4  0 
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APPENDIX  C 

STATISTICAL  ANALYSIS  OF  TOXICITY  TEST  DATA 
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***************£yj  EC50* ************************************* ******* 
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0  .450456D+03  .220000D+02  .700000D+03  .200000D+01 

1  .4044730+03  .2105050+02  .6664310+03  .1954070+01 

2  .4044490+03  .2108570+02  .6641760+03  .1940380+01 

3  .4044490+03  .2109000+02  .6639790+03  .1938760+01 

4  ,4044490+03  .2109040+02  .6639580+03  .1938600+01 

Oependent  variable  is  COCOONS 


Source  Sum-of 

-Squares 

df 

Mean-Square 

Regression 

7499.6 

3 

2499.9 

Residual 

404.4 

37 

10.9 

Total 

7904.0 

40 

Mean  corrected 

1999.1 

39 

Raw  R-square  (1-Residual/Total)  =  0.9 

Mean  corrected  R-square  { 1-Residual/Corrected)  =  0.8 

R (observed  vs  predicted)  square  =  0.8 

Wald  Confidence  Interval 
Lower  <  95%>  Upper 
18.5  23.7 

557.9  770.0 

1.1  2.8 


COCOONS 

COCOONS 

Case 

Observed 

Predicted 

Residual 

1 

27.0 

21.0 

6.0 

2 

18.0 

21.0 

-3.0 

3 

19.0 

21.0 

-2.0 

4 

20.0 

21.0 

-1.0 

5 

22.0 

20.3 

1.7 

6 

20.0 

20.3 

-0.3 

7 

22.0 

20.3 

1.7 

8 

23.0 

20.3 

2.7 

9 

20.0 

18.9 

1.1 

10 

11.0 

18.9 

-7.9 

11 

18.0 

18.9 

-0.9 

12 

19.0 

18.9 

0.1 

13 

16.0 

15.6 

0.4 

14 

17.0 

15.6 

1.4 

15 

13.0 

15.6 

-2.6 

16 

14.0 

15.6 

-1.6 

Parameter  Estimate  A.S.E.  Param/ASE 

G  21.1  1.3  16.2 

X  664.0  52.3  12.7 

B  1.9  0.4  4.7 
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17 

11.0 

12.4 

-1,4 

18 

15.0 

12.4 

2.6 

19 

10.0 

12.4 

-2.4 

20 

20.0 

12.4 

7.6 

21 

17.0 

8.9 

8.1 

22 

15.0 

8.9 

6.1 

23 

6.0 

8.9 

-2.9 

24 

11.0 

8.9 

2,1 

25 

9.0 

10.0 

-1,0 

26 

8.0 

10.0 

-2,0 

27 

9.0 

10.0 

-1.0 

28 

8.0 

10.0 

-2.0 

29 

5.0 

8.0 

-3.0 

30 

3.0 

8.0 

-5.0 

31 

5.0 

8.0 

-3.0 

32 

9.0 

8.0 

1.0 

33 

2.0 

4.6 

-2.6 

34 

5.0 

4.6 

0.4 

35 

5.0 

4.6 

0.4 

36 

8.0 

4.6 

3.4 

37 

1.0 

2.2 

-1.2 

38 

1.0 

2.2 

-1.2 

39 

3.0 

2.2 

0.8 

40 

1.0 

2.2 

-1.2 

Asymptotic  Correlation  Matrix 

G 


G 

X 

B 


1.0 

-0.8 

-0.6 


of  Parameters 
X 


1.0 

0.7 


B 


1.0 


Data,  estimates  and  residuals  have  been  saved. 


**************barium  ew  cocoon 


EC20* ****** ************* 


Iteration 
No.  Loss 

0  .111422D+04 

1  .568539D+03 

2  .461004D+03 

3  .4060070+03 

4  .4044540+03 

5  .4044490+03 

6  .4044490+03 

7  .4044490+03 


G 

.2200000+02 

.2227270+02 

.2043880+02 

.2111720+02 

.2110560+02 

.2109250+02 

.2109060+02 

.2109050+02 


X 

.2500000+03 

.2556240+03 

.4054960+03 

.3603370+03 

.3685680+03 

.3698840+03 

.3700030+03 

.3700150+03 


B 

.2000000+01 

.1117060+01 

.1836160+01 

.1917500+01 

.1929310+01 

.1937890+01 

.1938510+01 

.1938570+01 


Oependent  variable  is  COCOONS 


Source  Sum-of- 

•Squares 

df 

Mean-Square 

Regression 

7499.6 

3 

2499.9 

Residual 

404.4 

37 

10.9 

Total 

7904.0 

40 

Mean  corrected 

1999.1 

39 
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Raw  R-square  (l-Residual/Total) 

Mean  corrected  R-square  ( 1-Residual/Gorrected)  = 
R (observed  vs  predicted)  square  « 


Parameter 

G 

X 

B 


Wald  Confidence  Interval 
Estimate  A.S.E.  Param/ASE  Lower  <  95%>  Upper 

21-1  1.3  16.2  18.5  23.7 

370.0  69.1  5.4  230.1  509.9 

1-9  0.4  4.7  1.1  2.8 


1 

27.0 

21.0 

6.0 

2 

18,0 

21.0 

-3.0 

3 

19.0 

21.0 

-2.0 

4 

20.0 

21.0 

-1.0 

5 

22.0 

20.3 

1.7 

6 

20.0 

20.3 

-0.3 

7 

22.0 

20.3 

1.7 

8 

23.0 

20.3 

2.7 

9 

20.0 

18.9 

1.1 

10 

11.0 

18.9 

-7.9 

11 

18.0 

18.9 

-0.9 

12 

19.0 

18.9 

0.1 

13 

16.0 

15.6 

0.4 

14 

17.0 

15.6 

1.4 

15 

13,0 

15.6 

-2.6 

16 

14.0 

15.6 

-1.6 

17 

11.0 

12.4 

-1.4 

18 

15.0 

12.4 

2.6 

19 

10.0 

12.4 

-2.4 

20 

20.0 

12.4 

7.6 

21 

17.0 

8.9 

8.1 

22 

15.0 

8.9 

6.1 

23 

6.0 

8.9 

-2.9 

24 

11.0 

8.9 

2.1 

25 

.9.0 

10.0 

-1.0 

26 

8.0 

10.0 

-2.0 

27 

9.0 

10.0 

-1.0 

28 

8.0 

10.0 

-2.0 

29 

5.0 

8.0 

-3.0 

30 

3.0 

8.0 

-5.0 

31 

5.0 

8.0 

-3.0 

32 

9.0 

8.0 

1.0 

33 

2.0 

4.6 

-2.6 

34 

5.0 

4.6 

0.4 

35 

5.0 

4.6 

0.4 

36 

8.0 

4.6 

3.4 

37 

1.0 

2.2 

-1.2 

38 

1.0 

2.2 

-1.2 

39 

3.0 

2.2 

0.8 

40 

1.0 

2.2 

-1.2 

Asymptotic 

Correlation  Matrix 

of  Parameters 

G 

X 

B 

G 

1.0 

X 

-0.8 

1.0 

B 

-0.6 

0.9 

1.0 
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Data,  estimates  and  residuals  have  been  saved. 


GRAPH  MODEL  FOR  COCOONS: 

graph 

begin 

plot  juveniles*concentr  /  titles'*,  xlab=’RDX  mg/kg*,  ylab=*NUMBER  OF 
JUVENILES', 

xmax=150,  xmin=0,  ymax=30,  ymin=0 

fplot  y=17 . 732*exp  (  (log  ( .  5)  )  *  (concentr/2 . 389)  ^0 . 56) ;  xmin=0,  xmax==150,  xlab= 
ymin=0,  ylab='', 
ymax=30 

end 

SYSTAT  Rectangular  file  H:\SIMINI.SYD, 

created  Tue  Apr  09,  2002  at  10:40:18,  contains  variables: 

BA  COCOONS  ESTIMATE  RESIDUAL 

****************g.pg^  plot************************ 


Stem  and  Leaf  Plot  of  variable:  RESIDUAL,  N 


Minimum : 

-7,9 

Lower  hinge: 

-2.0 

Median: 

-0.9 

Upper  hinge: 

1.6 

Maximum : 

8.1 

40 


-7  8 

*  *  *  Outside  Values  *  *  * 

-4  9 

-3  0 

-2  H  998554000 
-1  H  54111000 
-0  M  82 
0  14448 

1  H  01477 

2  157 

3  4 

4 

5  9 

6  1 

*  *  *  Outside  Values  *  *  * 

7  5 

8  1 


ANOVA  for  Ba  effect  on  E.  fetlda  Cocoons 

TUE  4/09/02  2:48:26  PM 

SYSTAT  VERSION  7.0.1 
COPYRIGHT  (C)  1997,  SPSS  INC. 

Welcome  to  SYSTAT! 

SYSTAT  Rectangular  file  H:  \SIMINI-1\METPSSNL\BA\BADAT,SYD, 
created  Tue  Apr  09,  2002  at  11:05:26,  contains  variables: 
BA  COCOONS 
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Effects  coding  used  for  categorical  variables  in  model. 

Categorical  values  encountered  during  processing  are: 

BA  (10  levels) 

34,  153,  258,  433,  578,  689,  744, 

791,  1000,  1222 

Dep  Var:  COCOONS  N:  40  Multiple  R:  0.9  Squared  multiple  R:  0.9 

-1 

Estimates  of  effects  B  =  (X'X)  X'Y 


COCOONS 


CONSTANT 

12.2 

BA 

34 

8.9 

BA 

153 

9.6 

BA 

258 

4.9 

BA 

433 

2.9 

BA 

578 

1.9 

BA 

689 

"3.7 

BA 

744 

0.1 

BA 

791 

-6.7 

BA 

1000 

-7.2 

Analysis  of  Variance 


Source 

Sum-of -Squares 

df 

Mean-Square 

F-ratio 

BA 

1710.6 

9 

190.1 

19.8 

Error 

288.5 

30 

9.6 

Least  squares  means. 


LS  Mean 

SE 

N 

BA 

=34 

21.0 

1.6 

4 

BA 

=153 

21.8 

1.6 

4 

BA 

=258 

17.0 

1.6 

4 

BA 

=433 

15,0 

1.6 

4 

BA 

=578 

14.0 

1.6 

4 

BA 

=689 

8.5 

1.6 

4 

BA 

=744 

12.3 

1,6 

4 
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BA 

BA 

BA 


=791 

=1000 

=1222 


5.5 
5.0 

1.5 


1.6 

1.6 

1.6 


4 

4 

4 


Durbin-Watson  D  Statistic  2.127 

First  Order  Autocorrelation  -0.127 
COL/ 

ROW  BA 

1  34 

2  153 

3  258 

4  433 

5  578 

6  689 

7  744 

8  791 

9  1000 
10  1222 

Using  least  squares  means. 

Post  Hoc  test  of  COCOONS 


Using  model  MSE  of 

9.617  with  30  df. 

Matrix  of  pairwise 

mean  differences: 

1 

2 

1 

0.0 

2 

0.8 

0.0 

3 

-4.0 

-4.8 

4 

-6.0 

-6.8 

5 

-7.0 

-7.8 

6 

-12.5 

-13.3 

7 

-8.7 

-9.5 

8 

-15.5 

-16.3 

9 

-16.0 

-16.8 

10 

-19.5 

-20.3 

6 

7 

6 

0.0 

7 

3,8 

0.0 

8 

-3.0 

-6.8 

9 

-3.5 

-7.3 

10 

-7.0 

-10.8 

Fisher’s  Least-Significant-Difference  Test. 

Matrix  of  pairwise 

comparison  probabilities 

1 

2 

1 

1.0 

2 

0.7 

1.0 

3 

0.1 

0.0 

3  4  5 


0.0 

-2.0 

0.0 

-3.0 

-1.0 

0.0 

-8.5 

-6.5 

-5.5 

-4.8 

-2.8 

-1.8 

-11.5 

-9.5 

-8.5 

-12.0 

-10.0 

-9.0 

-15.5 

-13.5 

-12,5 

8 

9 

10 

0.0 

-0.5 

0.0 

-4.0 

-3.5 

o 

o 

3  4  5 

1.0 
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4  0.0  0.0 

5  0.0  0.0 

6  0.0  0.0 

7  0.0  0.0 

8  0.0  0.0 

9  0.0  0.0 

10  0.0  0.0 

6  7 

6  1.0 

7  0.1  1.0 

8  0.2  0.0 

9  0.1  0.0 

10  0.0  0.0 


0.4 

1.0 

o 

to 

0.7 

1.0 

o 

o 

0.0 

0.0 

o 

o 

0.2 

0.4 

o 

o 

0.0 

0.0 

o 

o 

0.0 

0.0 

o 

o 

0.0 

9 

0.0 

10 

1.0 

0.8 

1.0 

0.1 

0.1 

40  cases  have  been  saved  into  a  SYSTAT  file 


ANOVA  for  Ba  offoct  on  E.  fatida  Adults 


WED  4/10/02  12:19:54  PM 

SYSTAT  VERSION  7.0.1 
COPYRIGHT  (C)  1997,  SPSS  INC. 

Welcome  to  SYSTAT! 

SYSTAT  Rectangular  file  H:\SIMINI~l\ECOSSL\BABEMNSB\BA\BANAC.SYD, 
created  Thu  Aug  16,  2001  at  15:39:54,  contains  variables: 

CONC  SURVIVORS 

Effects  coding  used  for  categorical  variables  in  model. 

Categorical  values  encountered  during  processing  are: 

CONC  (9  levels) 

3^'  1000'  1124,  1222,  1348,  1585,  2000, 

2194,  2697 

Dep  Var:  SURVIVORS  N:  36  Multiple  R:  1.0  Squared  multiple  R:  0.9 


-1 


Estimates  of 

effects 

B  =  (X'X)  X'Y 

SURVIVORS 

CONSTANT 

3.2 

CONC 

34 

1.8 

CONC 

1000 

1.8 

CONC 

1124 

1.8 

CONC 

1222 

1.8 
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CONC 

1348 

1.8 

CONC 

1585 

0.6 

CONC 

2000 

-3.2 

CONC 

2194 

-3.2 

Source 

CONC 

Error 


Analysis  of  Variance 
Sum-of-Squares  df  Mean-Square 

188.9  8  23.6 

18.7  27  0.7 


F-ratio 

34.0 


0.0 


Least  squares  means. 


CONC 

=34 

CONC 

=1000 

CONC 

=1124 

CONC 

=1222 

CONC 

=1348 

CONC 

=1585 

CONC 

=2000 

CONC 

=2194 

CONC 

=2697 

LS  Mean  SE 

5.0  0.4 

5.0  0.4 

5.0  0.4 

5.0  0.4 

5.0  0.4 

3.7  0.4 

0.0  0.4 

0.0  0.4 

0.0  0.4 


N 


4 

4 

4 

4 

4 

4 

4 

4 

4 


Durbin-Watson  D  Statistic 
First  Order  Autocorrelation 
COL/ 

ROW  CONC 

1  34 

2  1000 

3  1124 

4  1222 

5  1348 

6  1585 

7  2000 

8  2194 

9  2697 

Using  least  squares  means. 
Post  Hoc  test  of  SURVIVORS 


2.167 

-0.083 
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Using  model  MSE  of 

0.694  with  27  df. 

Matrix  of  pairwise 

mean  differences: 

1  2 

1 

0.0 

2 

0.0  0.0 

3 

0.0  0.0 

4 

0.0  0.0 

5 

0.0  0.0 

6 

-1.3  -1.3 

7 

-5.0  -5.0 

8 

-5.0  -5.0 

9 

-5,0  -5.0 

6  7 

6 

0.0 

7 

-3.8  0.0 

8 

-3.7  0.0 

9 

-3.8  0.0 

Fisher’s  Least-Significant-Difference  Test. 

Matrix  of  pairwise 

comparison  probabilities 

1  2 

1 

1.0 

2 

1.0  1.0 

3 

1.0  1.0 

4 

1.0  1.0 

5 

1.0  1.0 

6 

0.0  0.0 

7 

0.0  0.0 

6 

0.0  0.0 

9 

0.0  0.0 

6  7 

6 

1,0 

7 

0.0  1.0 

8 

0.0  1.0 

9 

0.0  1.0 

,0.0 

0.0 

0.0 

-1.3 

-5.0 

-5.0 

-5.0 


0.0 

0.0 


0.0 

0.0 

-1.3 

-5.0 

-5.0 

-5.0 


0.0 


0.0 

-1.3 

-5.0 

-5.0 

-5.0 


1.0 

1.0 

1.0 

0.0 

0,0 

0.0 

0.0 


1.0 

1.0 


1.0 

1.0 

0.0 

0.0 

0.0 

0.0 


1.0 


1.0 

0.0 

0.0 

0.0 

0.0 


E.  fmtida  BERYLLIUM  NONLINEAR  REGRESSION  WITH  COCOONS 

TUE  4/09/02  1:59:40  PM 


SYSTAT  VERSION  7.0.1 
COPYRIGHT  (C)  1997,  SPSS  INC. 


****************EC50 


Welcome  to  SYSTAT! 

SYSTAT  Rectangular  file  H:\SIMINI~l\ECOSSL\BABEMNSB\BE\BECH2.SYD, 
created  Thu  Sep  21,  2000  at  16:08:20,  contains  variables: 

TRT  SURVIVAL  BIOMASS 
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Iteration 

No.  Loss  G  X  B 

0  .1416920+03  .8000000+01  .7000000+02  .2000000+01 

1  .1243450+03  .7821050+01  .5794980+02  .2313750+01 

2  .1188290+03  .7272510+01  .6438490+02  .3514490+01 

3  .1136270+03  .7324850+01  .6315280+02  .4968510+01 

4  .1133060+03  .7361020+01  .6312100+02  .5572380+01 

5  ,1133020+03  .7362130+01  .6306210+02  .5668370+01 

6  .1133020+03  .7361910+01  .6306600+02  .5671570+01 

7  .1133020+03  .7361900+01  .6306610+02  .5671700+01 


Oependent  variable  is  COCOONS 


Source  Sum- 

of-Squares 

df 

Mean-Square 

Regression 

745.7 

3 

248.6 

Residual 

113.3 

25 

4.5 

Total 

859.0 

28 

Mean  corrected 

403.0 

27 

Raw  R-square  (1-Residual/Total)  =  0.9 

Mean  corrected  R-square  (1-Residual/Corrected)  =  0.7 

R (observed  vs  predicted)  square  =  0.7 

Wald  Confidence  Interval 

Parameter  Estimate  A.S.E.  Param/ASE  Lower  <  95%>  Upper 

G  7.4  0.6  11.6  6.1  8.7 

X  63.1  6.0  10.6  50.8  75.4 

B  5.7  3.5  1.6  -1.6  13.0 


COCOONS  COCOONS 


Case 

Observed  * 

Predicted 

Residual 

1 

8,0 

7.4 

0.6 

2 

11.0 

7.4 

3.6 

3 

7.0 

7.4 

-0.4 

4 

3.0 

7.4 

-4.4 

5 

6.0 

7.3 

-1.3 

6 

7.0 

7.3 

-0,3 

.7 

6.0 

7.3 

-1.3 

8 

13.0 

7.3 

5.7 

9 

6.0 

7.3 

-1.3 

10 

10.0 

7.3 

2.7 

11 

6.0 

7.3 

-1.3 

12 

5.0 

7.3 

-2.3 

13 

5.0 

5.0 

0,0 

14 

5.0 

5.0 

0.0 

15 

8.0 

5.0 

3.0 

16 

2.0 

5,0 

-3.0 

17 

0.0 

0.3 

-0.3 

18 

0.0 

0.3 

-0.3 

19 

0.0 

0,3 

-0.3 

20 

1.0 

0.3 

0.7 

21 

1.0 

0.0 

1.0 

22 

0.0 

0.0 

0.0 

23 

3.0 

0.0 

3.0 

24 

0.0 

0.0 

0.0 

25 

0,0 

0.0 

0.0 
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26 

27 

28 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 


Asymptotic  Correlation  Matrix  of  Parameters 

G  X  B 

G  1.0 

X  -0.3  1,0 

B  -0.3  -0.1  1.0 

Data,  estimates  and  residuals  have  been  saved. 

******★***★★★*******★**★£020  COCOONS********************************* 
Iteration 

No.  Loss  G  X  B 

0  .190605D+03  .800000D+01  .500000D+02  .200000D+01 

1  .154972D+03  .8142460+01  .190620D+02  .1393610+01 

2  .1333700+03  .6953720+01  .4100540+02  .2378360+01 

3  .1193400+03  .7607410+01  .4269740+02  .3066020+01 

4  ,1163500+03  .7262520+01  .5212330+02  .4602740+01 

5  .1133890+03  .7372640+01  .5086370+02  .5236780+01 

6  .1133030+03  .7362280+01  .5161900+02  .5624600+01 

7  .1133020+03  .7362010+01  .5164000+02  .5671450+01 

8  .1133020+03  .7361900+01  .5164240+02  .5671650+01 

9  .1133020+03  .7361900+01  .5164250+02  .5671700+01 


Dependent  variable  is  COCOONS 


Source  Sum- 

of-Squares 

df 

Mean-Square 

Regression 

745.7 

3 

248.6 

Residual 

113.3 

25 

4.5 

Total 

859.0 

28 

Mean  corrected 

403.0 

27 

Raw  R-square  (1-Residual/Total)  =  0.9 

Mean  corrected  R-square  (1-Residual/Corrected)  =  0.7 

R (observed  vs  predicted)  square  =  0.7 


Parameter 

G 

X 

B 


Estimate 

7.4 

51.6 

5.7 


0.6 

7.5 

3.5 


Param/ASE 

11,6 

6.9 

1.6 


Wald  Confidence  Interval 
Lower  <  95%>  Upper 
6.1  8.7 

36,2  67.1 

-1.6  13.0 


COCOONS  COCOONS 

Case  Observed  Predicted  Residual 

1  8.0  7.4  0.6 

2  11.0  7.4  3.6 

3  7.0  7.4  -0.4 

4  3.0  7.4  -4.4 

5  6.0  7.3  -1.3 

6  7.0  7.3  -0.3 

7  6.0  7.3  -1.3 

8  13.0  7,3  5.7 

9  6.0  7.3  -1.3 
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10 

10.0 

7.3 

2.7 

11 

6.0 

7.3 

-1.3 

12 

5.0 

7.3 

-2.3 

13 

5.0 

5.0 

0.0 

14 

5.0 

5,0 

0.0 

15 

GO 

O 

5.0 

3.0 

16 

2.0 

5.0 

-3.0 

17 

o 

o 

0.3 

-0.3 

18 

o 

o 

0.3 

-0.3 

19 

0.0 

0.3 

-0.3 

20 

1.0 

0.3 

0.7 

21 

1.0 

0.0 

1.0 

22 

0.0 

0.0 

0.0 

23 

3.0 

0.0 

3.0 

24 

0.0 

0.0 

0.0 

25 

0.0 

0.0 

0.0 

26 

0.0 

0.0 

0.0 

27 

0.0 

0.0 

0.0 

28 

0.0 

0.0 

0.0 

Asymptotic  Correlation  Matrix  of  Parameters 

G  X  B 

G  1.0 

X  -0.4  1.0 

B  -0.3  0.8  1.0 

Data,  estimates  and  residuals  have  been  saved. 

SYSTAT  Rectangular  file  H:\simini.SYD, 

created  Tue  Apr  09,  2002  at  14:08:42,  contains  variables: 
BE  COCOONS  ESTIMATE  RESIDUAL 


SYSTAT  Rectangular  file  H:\simini.SYD, 

created  Tue  Apr  09,  2002  at  14:08:42,  contains  variables: 
BE  COCOONS  ESTIMATE  I^SIDUAL 

**************stem  and  leaf  plot************************** 

Stem  and  Leaf  Plot  of  variable:  RESIDUAL,  N  =  28 


Minimum: 

-4.4 

Lower  hinge: 

-0.8 

Median : 

0.0 

Upper  hinge: 

0.7 

Maximum : 

5.7 

-4  3 

*  *  *  Outside  Values  *  *  * 

-2  93 

-1  3333 

-0  M  33222000000 
0  M  00679 
1 

2  6 

*  *  *  Outside  Values  *  *  * 

2  9 

3  06 

5  6 


APPENDIX  C 


55 


ANOVA  £or  Be  effect  on  E,  fetida  Cocoons 

TUE  4/09/02  3:00:10  PM 

SYSTAT  VERSION  7.0.1 
COPYRIGHT  (C)  1997,  SPSS  INC. 

Welcome  to  SYSTAT! 

SYSTAT  Rectangular  file  H:\SIMINI~l\METANOVA\BA\BADAT.SYD, 
created  Tue  Apr  09,  2002  at  14:50:00,  contains  variables: 

BA  COCOONS 

40  cases  have  been  saved  into  a  SYSTAT  file 

Effects  coding  used  for  categorical  variables  in  model. 

Categorical  values  encountered  during  processing  are: 

BE  (7  levels) 

2-5,  24,  28,  57,  83,  121,  154 

Dep  Var:  COCOONS  N:  28  Multiple  R:  0.9  Squared  multiple  R:  0.7 

-1 

Estimates  of  effects  B  =  (X'X)  X'Y 


COCOONS 


CONSTANT 

4.0 

BE 

2.5 

3.2 

BE 

24 

4.0 

BE 

28 

2.7 

BE 

57 

1.0 

BE 

83 

-3.8 

BE 

121 

-3.0 

Analysis  of  Variance 


Source 

Sum-of-Squares 

df 

Mean-Square 

F-ratio 

BE 

296.7 

6 

49.5 

9.8 

Error 

106.3 

21 

5.1 
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***  WARNING  *** 

Case  8  is  an  outlier  (Studentized  Residual  =  3.0) 

Durbin-Watson  D  Statistic  2.262 

First  Order  Autocorrelation  -0.134 
COL/ 

ROW  BE 

1  2.5 

2  24 

3  28 

4  57 

5  83 

6  121 

7  154 

Using  least  squares  means. 

Post  Hoc  test  of  COCOONS 


Using  model  MSE  of  5.060  with  21  df. 
Matrix  of  pairwise  mean  differences; 


1  2 

1  0.0 

2  0.7  0.0 

3  -0.5  -1.2 

4  -2.2  -3.0 

5  -7.0  -7.7 

6  -6.2  -7.0 

7  -7.2  -8.0 

6  7 

6  0.0 

7  -1.0  0.0 

Fisher's  Least-Significant-Difference  Test. 
Matrix  of  pairwise  comparison  probabilities; 

1  2 

1  1.0 

2  0.6  1.0 

3  0.8  0.4 

4  0.2  0.1 


0.0 

-1.7 

-6.5 

-5.7 

-6.7 


1.0 

0.3 


4  5 


0.0 

-4.7  0.0 
-4.0  0.7 
-5.0  -0.3 


4  5 


1.0 
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5 

6 
7 


0.0  0.0 

0.0  0.0 

0.0  0.0 

6  7 

6  1.0 

7  0.5  1.0 


0.0  0.0  1.0 

0.0  0.0  0.6 

0.0  0.0  0.9 


ANOVA  for  Be  effect:  on  E.  fetida  Adults 

WED  4/24/02  3:17:22  PM 

SYSTAT  VERSION  7.0.1 
COPYRIGHT  (C)  1997,  SPSS  INC. 

Welcome  to  SYSTAT! 

SYSTAT  Rectangular  file  C:\SIMINI~l\ECOSSL\BABEMNSB\BE\BEAC2.SYD, 
created  Tue  Apr  09,  2002  at  13:52:48,  contains  variables: 

TRT  SURVIVAL  BIOMASS 

32  cases  have  been  saved  into  a  SYSTAT  file 
32  cases  have  been  saved  into  a  SYSTAT  file 

Effects  coding  used  for  categorical  variables  in  model. 

Categorical  values  encountered  during  processing  are: 

TRT  (8  levels) 

2.5,  79,  83,  110,  144,  191,  308, 

380 

Dep  Var:  SURVIVAL  N:  32  Multiple  R:  1.0  Squared  multiple  R:  1.0 


-1 

Estimates  of  effects  B  =  (X'X)  X'Y 


SURVIVAL 


CONSTANT 

3.0 

TRT 

2.5 

2.0 

TRT 

79 

2.0 

TRT 

83 

00 

I— 1 

TRT 

110 

o 

00 

TRT 

144 

1.5 

TRT 

191 

-2.0 

TRT 

308 

-3.0 
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Analysis  of  Variance 


Source 

TRT 

Error 


Sum-of-Squares  df  Mean-Square  F-ratio 

143.5  7  20,5  75.7 

6.5  24  0.3 


0.0 


Least  squares  means . 


LS  Mean 

SE 

N 

TRT 

-2.5 

5.0 

0,3 

4 

TRT 

=79 

5.0 

0.3 

4 

TRT 

=83 

00 

0.3 

4 

TRT 

=110 

00 

0.3 

4 

TRT 

=144 

4.5 

0.3 

4 

TRT 

=191 

1.0 

0.3 

4 

TRT 

=308 

0.0 

0.3 

4 

TRT 

=380 

0.0 

0.3 

4 

***  WARNING  *** 

Case  13  is  an  outlier 

Durbin-Watson  D  Statistic  2.135 

First  Order  Autocorrelation  -0.067 
COL/ 

ROW  TRT 

1  2,5 

2  79 

3  83 

4  110 

5  144 

6  191 

7  308 

8  380 

Using  least  squares  means. 

Post  Hoc  test  of  SURVIVAL 


(Student! zed  Residual  =  3.3) 


Using  model  MSE  of  0.271  with  24  df. 
Matrix  of  pairwise  mean  differences: 


1 

2 

3 

4 

5 

6 


1  2 

0.0 

0.0  0,0 

-0.2  -0.2 

-1.2  -1.2 

-0.5  -0.5 

-4.0  -4.0 


3 


0.0 

-1.0 

-0.2 

-3.7 


4  5 


0.0 

0.8  0.0 

-2.7  -3.5 
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7  -5.0 

“5.0 

-4.7 

-3.7 

-4.5 

8  -5.0 

“5.0 

“4.8 

“3.8 

-4.5 

6  7 

8 

6  0.0 

7  -1.0 

0.0 

8  -1.0 

0.0 

0.0 

Fisher's  Least-Significant-Difference 

Test. 

Matrix  of  pairwise  comparison  probabilities: 

1  2 

3 

4 

5 

1  1.0 

2  1.0 

1.0 

3  0.5 

0.5 

1.0 

4  0.0 

0.0 

0.0 

1.0 

5  0.2 

0.2 

0.5 

0.1 

1.0 

6  0.0 

0.0 

0.0 

0.0 

0.0 

7  0.0 

0.0 

0.0 

0,0 

0.0 

8  0.0 

0.0 

0.0 

0.0 

0.0 

6  7 

8 

6  1.0 

7  0.0 

1.0 

8  0.0 

1.0 

1.0 

E.  fBtida.  MANGANESE  NONLINEAR  RE(^SSION  WITH  COCOONS 


*************  ********£.050*  *********************  ******* 

Iteration 


No.  Loss 

0  .369828D+03 

1  .254551D+03 

2  .2480430+03 

3  .2478310+03 

4  .2478310+03 

5  .2478310+03 

6  .2478310+03 


G 

.7000000+01 

.8380340+01 

.8624140+01 

.8642340+01 

.8637960+01 

.8638380+01 

.8638350+01 


X 

.7000000+03 

.1003780+04 

.9155210+03 

.9260760+03 

.9266460+03 

.9266030+03 

.9266070+03 


B 

.2000000+01 

.2692000+01 

.2943780+01 

.2912890+01 

.2921610+01 

.2920840+01 

.2920910+01 


Oependent  variable  is  COCOONS 


Source  Sum- 

•of“Squares 

df 

Mean-Square 

Regression 

960.2 

3 

320.1 

Residual 

247.8 

24 

10.3 

Total 

1208.0 

27 

Mean  corrected 

440.0 

26 

Raw  R-square  (1-Residual/Total)  = 

Mean  corrected  R-square  (1-Residual/Corrected)  = 
R (observed  vs  predicted)  square  = 


0.8 

0.4 

0.4 
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Parameter 

G 

X 

B 


Estimate 

8.6 

926.6 

2.9 


A.S.E. 

1.5 

164.5 

1.9 


Wald  Confidence  Interval 
Param/ASE  Lower  <  95% >  Upper 

5.6  5.4  11.8 

5.6  587.1  1266.1 

1.5  -1.0  6.9 


Case 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


COCOONS 

Observed 

11.0 

5.0 

9.0 

14.0 

12.0 

3.0 

6.0 

9.0 

6.0 

12.0 

2.0 

6.0 

13.0 

3.0 

3.0 

3.0 

5.0 

1.0 

4.0 

2.0 

6.0 

2.0 

1.0 

3.0 

1.0 

1.0 

1.0 


COCOONS 

Predicted 

8.6 

8.6 

8.6 

8.2 

8.2 

8.2 

8.2 

7.6 

7.6 

7.6 

7.6 

6.4 

6.4 

6.4 

6.4 
3.0 
3.0 
3.0 
3.0 

2.7 

2.7 

2.7 

2.7 

1.7 

1.7 

1.7 

1.7 


Residual 

2.4 
-3.6 

0.4 

5.8 

3.8 
-5.2 
-2.2 

1.4 

-1.6 

4.4 
-5 . 6 
-0.4 

6.6 

-3.4 

-3.4 

-0.0 

2.0 

-2.0 

1.0 

-0.7 

3.3 
-0.7 
-1.7 

1.3 
-0.7 
-0.7 
-0.7 


Asymptotic  Correlation  Matrix  of  Parameters 

G  X 

G  1.0 


X 

B 


-0.8 

-0.7 


1.0 

0.8 


Data,  estimates  and  residuals  have  been  saved. 


B 


1.0 


Iteration 
No.  Loss 

0  .306381D+03 

1  .254893D+03 

2  .247900D+03 

3  .247831D+03 

4  .247831D+03 

5  .247831D+03 

6  .247831D+03 


G 

.700000D+01 

.841129D+01 

.862533D+01 

.863251D+01 

.863864D+01 

.863832D+01 

.863835D+01 


X 

.5000000+03 

,7089340+03 

.6184740+03 

.6298590+03 

.6285270+03 

,6286030+03 

.6285960+03 


B 

.2000000+01 

.3119000+01 

.2831850+01 

.2928900+01 

.2920300+01 

.2920970+01 

.2920900+01 
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Dependent  variable  is  COCOONS 


Source  Sum- 

of-Squares 

df 

Mean-Square 

Regression 

960.2 

3 

320.1 

Residual 

247.8 

24 

10.3 

Total 

1208.0 

27 

Mean  corrected  ’ 

440.0 

26 

Raw  R-square  (1-Residual/Total) 

Mean  corrected  R-square  { 1-Residual/Corrected) 
R (observed  vs  predicted)  square 


Parameter 

G 

X 

B 


Case 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


Estimate 

8,6 

628.6 

2,9 


COCOONS 


A.S.E. 

1.5 

255.2 

1.9 


Param/ASE 

5.6 

2.5 

1.5 


COCOONS 


11.0 

8.6 

2.4 

5.0 

8.6 

-3.6 

9.0 

8.6 

0.4 

14.0 

8.2 

5.8 

12.0 

8.2 

3.8 

3.0 

8.2 

-5.2 

6.0 

8.2 

-2.2 

9.0 

7.6 

1.4 

6.0 

7.6 

-1.6 

12.0 

7.6 

4.4 

2.0 

7.6 

-5.6 

6.0 

6.4 

-0.4 

13.0 

6.4 

6.6 

3.0 

6.4 

-3.4 

3.0 

6.4 

-3.4 

3.0 

3.0 

-0.0 

5.0 

3.0 

2.0 

1.0 

3.0 

-2.0 

4.0 

3.0 

1.0 

2.0 

2.7 

-0.7 

6.0 

2.7 

3.3 

2.0 

2.7 

-0.7 

1.0 

2.7 

-1.7 

3.0 

1.7 

1.3 

1.0 

1.7 

-0.7 

1.0 

1.7 

-0.7 

1.0 

1.7 

-0.7 

Wald  Confidence  Interval 
Lower  <  95%>  Upper 
5.4  11.8 

101.8  1155.3 

-1.0  6.9 


Asymptotic  Correlation  Matrix  of  Parameters 

G  X 

G  1.0 

X  -0.8  1.0 

B  -0.7  1.0 
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Data,  estimates  and  residuals  have  been  saved, 
created  Wed  Apr  17,  2002  at  10:18:34,  contains  variables: 
MN  COCOONS  ESTIMATE  RESIDUAL 

Stem  and  Leaf  Plot  of  variable:  RESIDUAL,  N  ==  27 


Minimum: 

-5.6 

Lower  hinge: 

-1.8 

Median: 

-0.7 

Upper  hinge: 

1.7 

Maximum: 

6.6 

-5  51 

-4 

-3  633 

-2  10 
-1  H  65 
-0  M  7776630 
0  39 

1  H  249 

2  3 

3  38 

4  4 

5  8 

6  6 


ANOVA  for  Ifo  effect  on  E.  fetlda  Cocoons 

Effects  coding  used  for  categorical  variables  in  model. 

Categorical  values  encountered  during  processing  are: 

MN  (7  levels) 

94,  386,  528,  697,  1067,  1111, 

Dep  Var:  COCOONS  N:  27  Multiple  R:  0.7  Squared  multiple  R: 

-1 

Estimates  of  effects  B  =  (X’X)  X'Y 


COCOONS 


CONSTANT 

5.4 

MN 

94 

2.9 

MN 

386 

3.3 

MN 

528 

00 

MN 

697 

0.8 

MN 

1067 

CM 

CM 

1 

MN 

1111 

-2.7 
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4 


Analysis  of  Variance 


Source 

MN 

Error 


Sum-of-Squares  df  Mean-Square  F-ratio 

194.6  6  32,4  2.6 

245.4  20  12.3 


0.0 


Least  squares  means. 


MN 

=94 

LS  Mean 

8.3 

SE 

2.0 

N 

3 

MN 

=386 

8.8 

1.8 

4 

MN 

=528 

7.3 

1.8 

4 

MN 

=697 

6.3 

GO 

f-H 

4 

MN 

=1067 

3.2 

1.8 

4 

MN 

=1111 

2.8 

1.8 

4 

MN 

=1236 

1.5 

1.8 

4 

Durbin-Watson  D  Statistic  2.421 
First  Order  Autocorrelation  -0.225 
COL/ 

ROW  MN 

1  94 

2  386 

3  528 

4  697 

5  1067 

6  1111 

7  1236 

Using  least  squares  means. 

Post  Hoc  test  of  COCOONS 


Using  model  MSE  of  12.271  with  20  df. 
Matrix  of  pairwise  mean  differences: 


1 

2 

3 

4 

5 

6 
7 

6 

7 


0.0 

0.4 

-1.1 

-2.1 

-5.1 

-5.6 

-6.8 

0.0 

-1.3 


2 


0.0 

-1.5 

-2.5 

-5.5 

-6.0 

-7.3 

7 


0.0 


3 


5 


0.0 

-1.0 

-4.0 

-4.5 

-5.8 


0.0 

-3.0 

-3.5 

-4.8 


0.0 

-0,5 

-1.8 
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Fisher's  Least-Significant-Difference  Test. 
Matrix  of  pairwise  comparison  probabilities: 


1 

2 

3 

4 

5 

6 
7 

6 

7 


1 

l.Q 

0.9 

0.7 

0.4 

0.1 

0.0 

0.0 

6 

1.0 

0.6 


1.0 

0.6 

0.3 

0.0 

0.0 

0.0 


1,0 


3  4  5 


1.0 

0.7 

1.0 

0.1 

0.2 

1.0 

0.1 

0.2 

0.8 

0.0 

0.1 

0.5 

27  cases  have  been  saved  into  a  SYSTAT  file 
***WARNING*** 

The  file  C:\SIMINI~l\METANOVA\MN\MNCOC5.SYD  was  read  for  processing,  and 
its  contents  have  been  replaced  by  saving 
the  processed  data  into  it. 

27  cases  have  been  saved  into  a  SYSTAT  file 


ANOVA  for  Mn  effect  on  E.  fetida  Adults 

WED  4/17/02  11:03:19  AM 

SYSTAT  VERSION  7.0.1 
COPYRIGHT  (C)  1997,  SPSS  INC. 

Welcome  to  SYSTAT! 

SYSTAT  Rectangular  file  C:\SIMINI~l\METPSSNL\MN\MNCOC5.SYD, 
created  Wed  Apr  17,  2002  at  10:27:36,  contains  variables: 
mn  COCOONS 

27  cases  have  been  saved  into  a  SYSTAT  file 
32  cases  have  been  saved  into  a  SYSTAT  file 

Effects  coding  used  for  categorical  variables  in  model. 

Categorical  values  encountered  during  processing  are: 

MN  (8  levels) 

94,  326,  449,  611,  767,  1111,  1444, 

2222 

Dep  Var:  ADULTS  N:  32  Multiple  R:  1.0  Squared  multiple  R:  1.0 


-1 

Estimates  of  effects  B  =  (X'X)  X'Y 
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ADULTS 


CONSTANT 

4.4 

MN 

94 

0.6 

MN 

326 

0.6 

MN 

449 

0.6 

MN 

611 

0.6 

MN 

767 

0.6 

MN 

1111 

0.3 

MN 

1444 

0.3 

Analysis  of  Variance 


Source  Sum-of-Squares  df  Mean-Square  F-ratio  P 


MN 

Error 

61.5  7 

2.3  24 

8.8 

0.1 

93.7  0.0 

Least 

squares  means. 

LS  Mean 

SE 

N 

MN 

=94 

5.0 

0.2 

4 

MN 

=326 

5.0 

0.2 

4 

MN 

=449 

5.0 

0.2 

4 

MN 

=611 

5.0 

0.2 

4 

MN 

=767 

5.0 

0.2 

4 

MN 

=1111 

4.7 

0.2 

4 

MN 

=1444 

4.7 

0.2 

4 

MN 

=2222 

0.7 

0.2 

4 

— 

***  WARNING  *** 

Case  22  is  an  outlier 

Case  25  is  an  outlier 

Case  32  is  an  outlier 

Durbin-Watson  D  Statistic  2.250 
First  Order  Autocorrelation  -0.250 
COL/ 

ROW  MN 

1  94 

2  326 


( Student i zed  Residual  = 
(Studentized  Residual  = 
(Studentized  Residual  = 


-3.4) 

-3.4) 

-3.4) 
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3  449 

4  611 

5  767 

6  1111 

7  1444 

8  2222 

Using  least  squares  means. 
Post  Hoc  test  of  ADULTS 


Using  model  MSB  of  0.094  with  24  df. 
Matrix  of  pairwise  mean  differences: 


1 

2 

3 

4 

5 

1  0.0 

2  0.0 

0.0 

3  0.0 

0.0 

0.0 

4  0.0 

0.0 

0.0 

0.0 

5  0.0 

0.0 

0.0 

0.0 

0.0 

6  -0.2 

-0.2 

-0.2 

-0,2 

-0.3 

7  -0.2 

-0.2 

-0.2 

-0.2 

-0.3 

8  -4.2 

-4.2 

-4.2 

-4.2 

-4.3 

6 

7 

8 

6  0.0 

7  0.0 

0.0 

00 

1 

o 

-4.0 

0.0 

Fisher's  Least-Significant-Difference  Test. 

Matrix  of  pairwise  comparison  probabilities: 

1 

2 

3 

4 

5 

1  1.0 

2  1.0 

1.0 

3  1.0 

1.0 

1.0 

4  1.0 

1.0 

1.0 

1.0 

5  1.0 

1.0 

1.0 

1.0 

1.0 

6  0.3 

0.3 

0.3 

0.3 

0.3 

7  0.3 

0.3 

0.3 

0.3 

0.3 

8  0.0 

0.0 

0.0 

0.0 

0.0 

6 

7 

8 

6  1.0 

7  1.0 

1.0 

8  0.0 

0.0 

1.0 

***WARNING*** 

The  file  C:\SIMINI-l\METANOVA\MN\MNADULT5.SYD  was  read  for  processing,  and 
its  contents  have  been  replaced  by  saving 
the  processed  data  into  it, 

32  cases  have  been  saved  into  a  SYSTAT  file 
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E.fetida  ANTIMONY  NONLINEAR  REGRESSION  WITH  COCOONS 


THU  8/16/01  9:51:32  AM 

SYSTAT  VERSION  7.0.1 
COPYRIGHT  (C)  1997,  SPSS  INC. 

Welcome  to  SYSTAT! 


SYSTAT  VERSION  7.0.1 
COPYRIGHT  (C)  1997,  SPSS  INC. 

nonlin 

print=long 

model  cocoons=g*exp(  (log(l-.5)  )  *  (conc/x)  ''b) 

save  c:\simini  data\ecossl\babemnsb\sbresidsb4  /  resid 

estimate/  start  =  9,  75,  2  iter=200 


Iteration 


No. 

Loss 

G 

X 

B 

0 

.577175D+02 

.9000000+01 

.750000D+02 

.200000D+01 

1 

.569808D+02 

.8825240+01 

.683194D+02 

.101910D+01 

2 

.490480D+02 

.8742410+01 

.704113D+02 

.1354380+01 

3 

.489823D+02 

.8742580+01 

.  697020D+02 

.1366960+01 

4 

.489819D+02 

.8744200+01 

.696049D+02 

.1361630+01 

5 

.489818D+02 

.8744040+01 

.  696221D+02 

.1362650+01 

6 

.489818D+02 

.8744070+01 

.  696187D+02 

.1362450+01 

7 

.489818D+02 

.8744060+01 

,6961940+02 

.1362490+01 

Dependent  variable  is  COCOONS 

Source  Sum- 

of-Squares 

df  Mean-; 

Square 

Regression 

479.018 

3  159.673 

Residual 

48.982 

28 

1.749 

Total 

528.000 

31 

Mean 

corrected 

266.710 

30 

Raw  R-square  (1-Residual/Total)  =  0.907 

Mean  corrected  R-square  { 1-Residual/Corrected)  =  0.816 

R(observed  vs  predicted)  square  =  0.817 

Wald  Confidence  Interval 
Estimate  A.S.E.  Param/ASE  Lower  <  95% >  Upper 

8.744  0.659  13.264  7.394  10.094 

69.619  9.892  7.038  49.356  89.882 

1.362  0.316  4.306  0.714  2.011 

COCOONS  COCOONS 

Case  Observed  Predicted  Residual 

1  9.000  8.744  0.256 

2  8.000  8.744  -0.744 


Parameter 

G 

X 

B 
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3 

11.000 

8.744 

2.256 

4 

7.000 

8.744 

-1.744 

5 

5.000 

4.965 

0.035 

6 

6.000 

4.965 

1.035 

7 

5.000 

4.965 

0.035 

8 

4.000 

3.469 

0.531 

9 

4.000 

3.469 

0.531 

10 

1.000 

3.469 

-2.469 

11 

1.000 

3.469 

-2.469 

12 

0.0 

2.640 

-2.640 

13 

5.000 

2.640 

2.360 

14 

4.000 

2.640 

1.360 

15 

2.000 

2.640 

-0.640 

16 

3.000 

1.909 

1.091 

17 

2.000 

1.909 

0.091 

18 

3.000 

1.909 

1.091 

19 

4.000 

1.909 

2.091 

20 

2.000 

1.238 

0.762 

21 

0.0 

1.238 

-1.238 

22 

1.000 

1.238 

-0.238 

23 

2.000 

1.238 

0.762 

24 

0,0 

0.711 

-0,711 

25 

.  0.0 

0.711 

-0.711 

26 

0.0 

0.711 

-0.711 

27 

1.000 

0.711 

0.289 

28 

0.0 

0.349 

-0.349 

29 

0.0 

0.349 

-0.349 

30 

0.0 

0.349 

-0.349 

31 

0.0 

0,349 

-0.349 

Asymptotic  Correlation  Matrix  of  Parameters 

G  X 

G  1.000 

X  -0.547  1.000 

B  -0.288  0.793 


B 

1.000 


Data,  estimates  and  residuals  have  been  saved. 


Iteration 
No.  Loss 

0  .234717D+03 

1  .578881D+02 

2  .501699D+02 

3  .489842D+02 

4  .489819D+02 

5  .489818D+02 

6  .489818D+02 

7  .489818D+02 

8  .489818D+02 


G  X 

.900000D+01  .750000D+02 
.846243D+01  .194391D+02 
.874431D+01  .320674D+02 
.873859D+01  .305325D+02 
.874440D+01  .302514D+02 
.874400D+01  .303102D+02 
.874408D+01  .302985D+02 
.874406D+01  .303007D+02 
.874406D+01  .303003D+02 


B 

.200000D+01 

.934868D+00 

.136049D+01 

.137175D+01 

.136057D+01 

.136286D+01 

.136241D+01 

.136250D+01 

.136248D+01 


Dependent  variable  is  COCOONS 


Source 

Sum-of-Squares 

df 

Mean-Square 

Regression 

479.018 

3 

159.673 

Residual 

48.982 

28 

1.749 

Total 

528.000 

31 
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Mean  corrected 


266,710 


30 


Raw  R-square  (1-Residual/Total) 

Mean  corrected  R-square  { 1-Residual/Corrected)  = 
R (observed  vs  predicted)  square 


0.907 

0.816 

0.817 


Parameter 

G 

X 

B 


Estimate 

8.744 

30.300 

1.362 


Wald  Confidence  Interval 
A.S.E.  Param/ASE  Lower  <  95%>  Upper 

0-659  13.265  7.394  10.094 

9-633  3.145  10.568  50.033 

0-316  4.306  0.714  2.011 


COCOONS 


COCOONS 


1 

9.000 

8.744 

0.256 

2 

8.000 

8.744 

-0.744 

3 

11.000 

8.744 

2.256 

4 

7.000 

8.744 

-1.744 

5 

5.000 

4.965 

0.035 

6 

6.000 

4.965 

1.035 

7 

5.000 

4.965 

0.035 

8 

4.000 

3.469 

0.531 

9 

4.000 

3.469 

0.531 

10 

1.000 

3.469 

-2.469 

11 

1.000 

3.469 

-2.469 

12 

0.0 

2.640 

-2.640 

13 

5.000 

2.640 

2.360 

14 

4.000 

2.640 

1.360 

15 

2.000 

2.640 

-0.640 

16 

3.000 

1.909 

1.091 

17 

2.000 

1.909 

0.091 

18 

3.000 

1.909 

1.091 

19 

4.000 

1.909 

2.091 

20 

2.000 

1.238 

0.762 

21 

0.0 

1.238 

-1.238 

22 

1.000 

1.238 

-0.238 

23 

2.000 

1.238 

0.762 

24 

0.0 

0.711 

-0.711 

25 

0.0 

0.711 

-0.711 

26 

0.0 

0.711 

-0.711 

27 

1.000 

0.711 

0.289 

28 

0.0 

0.349 

-0.349 

29 

0.0 

0.349 

-0.349 

30 

.0.0 

0.349 

-0.349 

31 

0.0 

0.349 

-0.349 

Asymptotic  Correlation  Matrix  of  Parameters 

G  X 

G  1.000 

X  -0.420  1.000 

B  -0.288  0.962 


1.000 


Data,  estimates  and  residuals  have  been  saved. 
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>graph 


>begin 

>plot  cocoons*conc  /  title=**,  xlab='Sb  S04  concentration',  ylab=' NUMBER  OF 
COCOONS ' , 

>  xmax=250,  xmin=0,  ymax=12,  ymin=0 

> 

>fplot  y=8 .74*exp( {log( .5) ) * (conc/69. 62) ^1,36) ;  xmin^O,  xinax=250,  xlab='* 
yniin=0,  ylab='*,  yinax=12 


>enci 

>graph 

>begin 

>plot  cocoons*conc  /  title='*/  xlab='Sb  concentration',  ylab='NUMBER  OF 
COCOONS ' , 

>  xmax=250,  xmin=0,  yinax=12,  ymin=0 

> 

>graph 

>use  d:\nonlin\residsb3  /  resid 

SYSTAT  Rectangular  file  d: \\nonlin\residsb3 . SYD, 
created  Thu  Aug  16,  2001  at  10:06:06,  contains  variables: 

COCOONS  CONC  ESTIMATE  RESIDUAL 

>plot  residual*conc 

>plot  residual*estiinate 

> 

ANOVA  for  Sb  effect  on  E.  fetlda  Cocoons 

MON  9/18/00  2:34:10  PM 

SYSTAT  VERSION  7.0.1 
>ANOVA 

>CATEG0RY  TRT/MISS 
>COVAR 

>DEPEND  COCOONS  /  BONF 
>ESTIMATE 

Effects  coding  used  for  categorical  variables  in  model. 
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149 


179 


Categorical  values  encountered  during  processing  are: 

TRT  {8  levels) 

0,  60,  86,  104,  124, 

215 

Dep  Var:  COCOONS  N:  35  Multiple  R:  0.83  Squared  multiple  R:  0.69 

Analysis  of  Variance 

Source  Sum-of-Squares  df  Mean-Square  F-ratio  P 

TRT  230.81  7  32.97  8.66  0.00 

Error  102.79  27  3.81 


Durbin-Watson  D  Statistic  1.853 
First  Order  Autocorrelation  0.000 
COL/ 

ROW  TRT 
1  0 
2  60 

3  86 

4  104 

5  124 

6  149 

7  179 

8  215 

Using  least  squares  means. 

Post  Hoc  test  of  COCOONS 


Using  model  MSE  of  3.807  with  27  df. 
Matrix  of  pairwise  mean  differences: 


1 

2 

3 

4 

5 

1 

o 

o 

2 

-1.79 

o 

o 

3 

-4.63 

-2.83 

0.0 

4 

-4.38 

-2.58 

0.25 

o 

o 

5 

-4.13 

-2.33 

0.50 

0.25 

o 

o 

6 

CO 

00 

IT) 

1 

-4.08 

-1.25 

-1.50 

-1.75 

7 

00 

00 

-5.08 

-2.25 

-2.50 

-2.75 

8 

-7.13 

-5.33 

-2.50 

-2.75 

o 

o 

ro 

1 

6 

7 

8 

6 

0.0 

7 

-1.00 

0.0 

8 

-1.25 

-0.25 

o 

o 
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Bonferroni  Adjustment. 

Matrix  of  pairwise  comparison  probabilities: 


1 

2 

3 

4 

5 

1 

1.00 

2 

1.00 

1.00 

3 

0.02 

1.00 

1.00 

4 

0.03 

1.00 

1.00 

1.00 

5 

0.05 

1.00 

1.00 

1.00 

1.00 

6 

0.00 

0.30 

1.00 

1.00 

1.00 

7 

0.00 

0.06 

1.00 

1.00 

1.00 

8 

0.00 

0.04 

1.00 

1.00 

1.00 

6 

7 

8 

6 

1.00 

7 

1.00 

1.00 

8 

1.00 

1.00 

1.00 

ANOVA  for  Sb  effect  on  E.  fetlda  Adults 

WED  9/20/00  3:04:44  PM 

SYSTAT  VERSION  7.0.1 
COPYRIGHT  (C)  1997,  SPSS  INC. 

>ANOVA 

>CATEGORY  SBCONC 
>COVAR 

>DEPEND  SURVIVORS  /  BONF 
>ESTIMATE 

Effects  coding  used  for  categorical  variables  in  model. 

Categorical  values  encountered  during  processing  are: 

SBCONC  (7  levels) 

0,  245,  318,  412,  537,  617,  697 

Dep  Var:  SURVIVORS  N:  32  Multiple  R:  0.68  Squared  multiple  R:  0.46 

Analysis  of  Variance 


Source 

Sum-of  “Squares 

df 

Mean-Square 

F-ratio 

P 

SBCONC 

9.00 

6 

1.50 

3.57 

0.01 

Error 

10.50 

25 

0.42 
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***  WARNING  *** 

Case  30  is  an  outlier 

Case  32  is  an  outlier 

Durbin-Watson  D  Statistic  2.036 
First  Order  Autocorrelation  -0.125 
COL/ 

ROW  SBCONC 
1  0 

2  245 

3  318 

4  412 

5  537 

6  617 

7  697 

Using  least  squares  means. 

Post  Hoc  test  of  SURVIVORS 


(Studentized  Residual 
{Studentized  Residual 


3.10) 

-3.10) 


Using  model  MSE  of  0.420  with  25  df. 
Matrix  of  pairwise  mean  differences: 


1 

2 

3 

4 

5 

1 

0.0 

2 

0.00 

0.0 

3 

0.0 

0.00 

0.0 

4 

-0.25 

-0.25 

-0.25 

0.0 

5 

-0.25 

-0.25 

-0.25 

0.0 

0.0 

6 

-1.00 

-1.00 

-1.00 

-0.75 

-0.75 

7 

-1.50 

-1.50 

-1.50 

-1.25 

-1.25 

6 

7 

6 

o 

o 

7 

o 

IT) 

O 

o 

o 

Bonferroni  Adjustment. 

Matrix  of  pairwise  comparison  probabilities: 


1  2  3  4  5 


1 

1.00 

2 

1.00 

1,00 

3 

1.00 

1.00 

1.00 

4 

1.00 

1.00 

1.00 

1.00 

5 

1,00 

1.00 

1.00 

1.00 

1.00 

6 

0.39 

0.81 

0.81 

1.00 

1.00 

7 

0.02 

0.07 

0.07 

0.24 

0.24 

6 

7 

6 

1.00 

7 

1.00 

1.00 

APPENDIX  C 


74 


